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1. 


INTRODUCTION 


Vibration  is  a  phenomenon  specially  appealing  to  a 
number  of  engineers  due  partly  to  a  clearly  defined,  exten¬ 
sive  and  intriguing  body  of  underlying  mathematical  analysis. 
From  the  late  1800's  when  attention  was  first  paid  to  the 
problem  of  ship  hull  vibration  this  body  of  analysis  has  been 
expanded  by  a  large  number  of  active  researchers  to  provide 
the  current  comprehensive  level  of  understanding  of  the  subject. 

It  is  apparent  that  vibration  of  the  ship  hull  is  a 
source  of  concern  for  the  naval  architect  (and  shipowner) . 

A  passenger  vessel  with  vibration  amplitudes  above  a  certain 

\ 

threshold  would  have  few  repeat  passengers.  Fatigue  of 
structural  components  due  to  the  cyclic  stresses  is  a  factor 
to  be  considered.  The  conventional  method  of  determining  the 
strength  of  the  hull  girder  does  not  explicitly  take  into 
account  the  cyclic  stress  associated  with  hull  vibration. 

It  would  seem  necessary  to  determine  what  these  stresses 
might  be,  how  they  vary  with  time  and  the  effect  of  their 
addition  to  stresses  already  present  through  ship  loading 
and  the  quasistatic  effect  of  wave  support.  The  sophisticated 
equipment  on  board  warships  experience  problems  in  an  environ¬ 
ment  of  excessive  vibration.  It  is  interesting  to  note  that 
with  continued  refinement  of  such  equipment  as  navigation 
instruments,  radar,  sonar,  fire  control  computers,  communica¬ 
tions  gear  and  electronic  warfare  devices  the  requirement  for 
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vibration  free  foundations  becomes  more  stringent.  There 
are  dangers  to  shipping  from  acoustic  and  pressure  operated 
offensive  weapons.  Attention  has  been  focused  on  noise 
emitted  by  hulls  and  hull  vibration  in  general.  The  "acoustic 
signatures"  of  vessels  play  an  important  part  in  tracking 
foreign  vessels, both  surface  and  subsurface,  by  passive  means. 
It  appears  that  the  problem  of  vibration  control  is  becoming 
more  important. 

As  stated,  the  body  of  analysis  is  large  but  it  is 
apparent  that  the  complexity  of  hull  structures  in  both  form, 
material  and  construction  and  the  complexity  of  three-dimen¬ 
sional  hydrodynamics  combine  to  limit  the  application  of  the 
analysis  to  reality.  It  is  not  surprising  then  that  the  pro¬ 
blems  of  hull  vibration  still  remain  and  that  work  continues 
on  new  methods  of  mathematical  treatment,  on  models  and  on 
full  size  ships  where  the  adequacy  of  the  technical  progress 
can  be  shown. 

This  thesis  is  concerned  with  one  of  the  two  broad 
classes  of  ship  hull  vibration,  that  is  synchronous  or 
resonant  vibration  as  opposed  to  local  vibration.  In  this 
case  the  entire  hull  acts  like  a  beam  vibrating  under  some 
excitation  from  machinery,  the  propellor  or  waves.  In 
serious  cases  the  motion  can  clearly  be  seen  by  sighting 
along  the  length  of  the  ship.  The  importance  of  a  resonant 
vibration  is  dependent  on  the  conditions  of  its  occurence. 
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If  the  engine  at  normal  operating  speed  excites  the  hull  then 
a  serious  problem  exists.  If,  however,  this  resonance  occurs 
at  an  engine  speed  that  will  be  only  infrequently  maintained, 

v 

then  there  is  a  problem  of  considerably  less  significance. 

Similarly,  if  a  vessel  is  expected  to  operate  in  certain 
known  sea  states  then  it  is  necessary  to  avoid  an  excessive 
resonant  condition  by  proper  hull  design. 

The  other  class  of  hull  vibration  is  a  local  effect. 

Here  specific  parts  of  the  ship  or  certain  fittings  such  as 
masts,  stacks,  superstructure  or  a  panel  of  plating  are  set 
into  motion.  Depending  on  the  point  of  view  local  vibration 
can  also  include  the  undesired  motion  of  items  as  small  as 
light  fixtures  or  stanchion  chains .  Local  vibration  can  have 
the  same  effects  as  hull  resonance,  that  is  personal  discom¬ 
fort  and  structural  and  equipment  failures.  This  type  of 
vibration  is  not  as  readily  predicted  and  analyzed  as  whole 
hull  vibration  and  is  treated  generally  after  the  fact  by 
locating  and  removing  the  excitation  or  adding  local 
stiffening,  pillars  and  similar  structures.  Obviously  in 
many  cases  the  excitation  cannot  be  eliminated  and  additional 
structure  is  -not  feasible.  Then  the  solution  is  the  removal 
of  the  object  or  relaxation  of  the  requirements. 

It  is’  apparent  that  the  hull  of  a  ship  will  vibrate 
only  if  an  external  force  is  applied  to  it.  When  this  force 
is  from  slamming,  the  resulting  shudder  can  be  picked  up  by 
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accelerometers  and  can  be  shown  to  be  transient  in  nature. 

A  similar  force  can  be  generated  by  rapidly  releasing  a 
haused  anchor  or  rapidly  braking  a  free  falling  one.  The 
transient  response  and  information  which  can  be  obtained  from 
it  will  be  discussed  later.  Other  sources  of  excitation  can 
produce  continuous  vibration  which  is  somewhat  more  insidious 
due  to  its  gradual  deleterious  action.  Such  sources  include 
those  whose  origin  is  some  out-of-balance  periodic  force  in 
the  diesel  main  engine,  auxiliary  machinery,  shafting  or 
propellers. 

It  can  therefore  be  reduced  in  any  given  instance  by 
proper  attention  to  balancing  during  construction.  Diesel 
engine  balancing  is  a  well  developed  art  which  depends  on  the 
type  of  engine,  number  of  cylinders,  auxiliaries  which  may  be 
run  from  the  main  engine  and  many  other  factors.  Properly 
operating  turbine  equipment  runs  at  angular  velocities  well 
above  the  frequencies  of  concern  here. 

The  size  of  shafting  and  propellers  makes  perfect 
balancing  impossible.  In  addition  other  once  per  revolution 
forces  can  be  created  if  one  blade  of  the  propeller  has  a 
different  shape  than  the  other  blades.  Even  with  a  perfect 
propeller,  the  blades  work  in  a  mixed  wake  behind  the  hull  and 
the  force  on  each  varies  continuously  throughout  each  revolu¬ 
tion.  These  varying  forces  are  conveniently  divided  into  two 
types  according  to  the  method  of  transmission  to  the  hull. 


"Propeller"  forces  are  transmitted  through  the  shafting  to 
struts  and  bearings.  "Surface"  forces  are  transmitted 
through  the  water  to  the  stern.  Both  types  have  a  frequency 
of  propeller  angular  velocity  times  the  number  of  blades  and 
are  referred  to  as  blade-rate  forces.  Extensive  work  has  been 
carried  out  and  resultantly  the  magnitude  of  these  forces  can 
be  calculated  quite  accurately.  For  example,  the  propeller 
forces  can  be  determined  from  lifting  surface  theory  assuming 
the  hull  induced  velocity  field  in  the  propeller  plane  is 
known.  With  the  requirement  that  the  fluid  particles  should 
remain  in  contact  with  the  propeller  blades,  it  follows 
that  the  velocity  components  normal  to  the  surface  of  the 
blades  due  to  the  hull  induced  wake  should  be  equal  to  the 
corresponding  components  of  the  propeller  induced  fluid 
velocities.  This  requirement  can  be  expressed  as  a  three- 
dimensional  integral  equation.  The  numerical  computation  is 
formidable  but  computer  solutions  yield  the  pressure  distri¬ 
bution  and  the  blade  lift.  Integration  yields  the  forces. 

It  is  interesting  to  note  that  due  to  the  dynamic  behavior 
of  the  shaft,  bearing  reactions  can  be  an  order  of  magnitude 
higher  than  the  forces  generated  at  the  propeller.  All  of 
the  forces  can  be  kept  within  tolerable  limits  by  careful 
manufacture  and  proper  attention  to  stern  and  appendage 
design  and  particularly  to  blade  tip  clearance.  Some 
hydrodynamic  disturbing  forces  have  to  be  accepted,  however, 
as  part  of  the  price  of  mechanical  propulsion. 
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The  excitation  due  to  wave  passage  is  at  present  an 


active  area  of  research.  The  circular  frequency  of  encounter 
of  waves,  m  ,  is  dependent  on  the  wave  spectrum  of  the  sea 
surface  and  the  speed  of  the  vessel.  Even  for  the  swiftest 
of  destroyers  the  amount  of  energy  available  from  excitation 
of  even  the  lowest  natural  frequency  is  very  small,  the 
preponderance  of  wave  energy  being  at  lower  frequencies . 

Ships  of  more  recent  vintage,  however,  are  longer  with  some 
lengths  over  1,000  feet  and  natural  frequencies  for  the  two 
node  mode  of  between  3  and  5  radians/sec.  There  is  ample 
energy  for  excitation  in  this  range  and  numerous  studies  have 
concluded  that  ''springing”  as  two  node  vibration  is  called 
is  a  significant  addition  to  midship  stress  levels.  [1]  [12] 
[16] [30] 

When  engines  were  first  put  on  ships  hull  vibration 
problems  arose.  Shortly  thereafter  an  effort  was  made  to 
predict  the  frequency  for  the  two  node  mode  of  vibration. 

As  a  result  of  very  early  work  by  Otto  Schlich  and  by  F.W. 
Todd,  empirical  formulae  were  developed  which  closely 
correlated  with  observed  frequencies.  Schlich  proposed  a 
modification  to  the  natural  frequency  formula  for  a  uniform 
free-free  beam  which  is 
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where  32L  =  4,730,  E  is  the  Young's  modulus,  I  is  the  section 
moment  of  inertia,  m  is  the  beam  mass  per  unit  length,  and 
L  is  the  length.  This  formula  will  be  explained  more  fully 
later.  Then 


to,  -  N  (2) 

2  /  AL3 

4 

where  the  units  are  tons  for  A,  feet  for  L,  ft  for  I,  and 
radians/sec  for  co2.  Schlick  gave  the  following  values  for  N: 


ships  with  fine  lines . 1369. 

passenger  liners . 1252. 

cargo  ships . 1116. 


Another  early  effort  which  was  twin  steps  in  the 
directions  of  simplicity  and  accuracy  by  F.W.  Todd  was  based 
on 

I  =  CBD3  (3) 


where  a  is  the  beam,  D  is  the  depth,  and  C  a  coefficient 
dependent  on  L/D  and  the  ship  type.  In  the  Todd  formula, 
the  virtual  mass  was  used  which  was  a  function  of  the  beam 
to  draft  ratio  B/T.  So 


0) 


2 


(4) 
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where 

d.  =  A (1. 2  +  B/3T)  (5) 

Values  for  3  are  listed  in  [32] . 

In  the  treatment  which  follows  the  natural  frequencies 
of  a  vessel  are  determined  from  offset  data  and  other  infor¬ 
mation.  The  process  is  an  iterative  one  and  can  be  expedited 
by  an  appropriate  choice  of  a  starting  frequency.  The  choice 
of  this  value  is  from  a  formula  devised  by  L.C.  Burrill. 

w,  -  80,  OOOir  [ — * - - - ]1/2  (6) 

^  ALJ(1+B/2T) (1+r  ) 

S 

where 

_  3.5D2[3(B/D)3  +  9(B/D)2+6(B/D)+1.2] 
r  - 

S  u  ( 3B/D+1) 

The  Taylor  shear  correction  term  rg  is  dimensionless. 

The  analysis  of  hull  vibration  has  proceeded  far  beyond 
that  which  produced  the  preceding  empirical  formulations  and 
this  work  attempts  to  explore  the  advances  in  considerations 
of  damping,  excitation,  added  mass  and  the  use  of  the  computer. 
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TECHNICAL  TREATMENT 


A.  Uniform  Beam  Vibration 

As  previously  stated,  the  mathematical  base  for  the 
analysis  of  vibrating  bodies  is  large  and  acceptance  of 
certain  broad  assumptions  allows  complex  problems  to  be 
idealized  and  working  solutions  obtained.  As  more  realistic 
assumptions  are  made  more  intimate  knowledge  of  the  process 
is  obtained.  The  price,  of  course,  is  vastly  more  complex 
solution  techniques. 

Consider  a  beam  of  length  L.  The  vertical  displacement 
at  any  time,  t,  is  y(x,t)  where  x  is  the  longitudinal  axis 
of  the  beam.  The  total  vertical  force  per  unit  length  is 
f(x,t).  The  system  parameters  are  the  mass  per  unit  length 
m(x)  and  the  flexural  rigidity  EI(x)  where  E  is  Young's 
modulus  of  elasticity  and  l(x)  is  the  cross  section  area 
moment  of  inertia  about  the  neutral  axis  of  bending  at  any 
position  x  along  the  beam.  Consider  an  element  of  length 
dx  of  the  beam.  V(x,t)  and  M(x,t)  are  the  shearing  force 
and  bending  moment.  The  simplifications  to  be  included  at 
this  point  are  1)  shear  deformation  is  small  in  relation  to 
bending  deformation  and  2)  rotary  inertia  effects  are  small. 

These  simplifications  are  valid  if  the  beam's  length 
to  depth  ratio  L/D  is  about  10  or  greater  and  if  deflections 
and  slopes  are  small.  Most  ships  underway  fall  into  this 
category.  [14]  Another  major  assumption  is  that  damping  is 


small  enough  to  be  neglected.  This  is  sufficiently  accurate 
for  vibrating  ships  particularly  in  the  free  vibration 
analysis  which  follows.  [12]  The  small  damping  present 
serves  to  limit  vibration  amplitude  rather  than  modify 
frequency  or  vibration  pattern.  The  force  equation  of  motion 
in  the  vertical  direction  is  then 

2 

[V  (x, 1)  dx]  -  V  (x ,  t)  +f  (x / 1)  dx=m  (x)  dx3--X  — 

dx  3  c 

(1) 

The  moment  equation  of  motion,  ignoring  the  inertia  torque 
associated  with  the  element  rotation  is 

[M(x, t)  +  3M^ <*x]-M(x,t)  +  [V(x,t)+3-^x,t)  dx] dx 


+f  (x, t) dx  ^  = 


(2) 


Here  we  assumed  that  the  elemental  length  dx  is  small  enough 
so  that  f(x,t)  is  constant  along  the  element.  Equation  (2) 
reduces  to 

— +  V(x,t)  =  0  (3) 


Using  (3)  in  (1)  and  simplifying  results  in 


From  basic  strength  of  materials  it  is  known  that 


M(x,t)  =  EI(x)  9 

3x^ 


(5) 


(5)  into  (4)  yields 

~3-L[EI(x)  J+f  (X/t)  =  m(x) 

3x  3x  3t 


(6) 


which  is  the  differential  equation  for  the  flexural  vibration 
of  a  beam.  [25]  This  is  a  fourth  order  linear  nonhomogenous 
partial  differential  equation.  General  solution  theory  re¬ 
quires  a  number  of  boundary  conditions  equal  to  the  order  of 
the  equation.  Since  this  equation  must  be  integrated  four 
times  to  get  a  solution,  four  arbitary  constants  are  intro¬ 
duced,  thus  the  requirement  for  four  boundary  conditions. 

In  the  case  of  interest  the  beam  ends  are  free.  Thus  the 
shear  and  moment  at  the  ends  equal  zero.  From  equations  (3) 
and  (5) 


M(x,t)  =  EI(x)  -  =  o 


3x 


V(x,t)  =  -i-[EI(x)  3.  =  0 

3x  3x^ 


(7) 

(8) 


x  *  0,  L 
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The  next  analytical  step  is  consideration  of  free  vibration. 
The  equation  of  notion  reduces  to 


^jtEI(x) 

3xz  3x^ 


m(x) 


9zv(x,t) 

3t2 


(9) 


with  the  boundary  conditions  (7)  and  (8) .  This  forms  a 
classic  "boundary  value”  problem.  It  is  opportune  at  this 
point  to  explore  the  possibility  of  synchronous  motion.  This 
is  a  motion  in  which  the  general  shape  of  the  beam  displace¬ 
ment  does  not  change  with  time.  That  is  to  say,  every  point 
of  the  beam  executes  the  same  motion,  passing  through  the 
equilibrium  position  and  the  point  of  maximum  excursion  at 
the  same  time.  This  implies  that  the  deformation  g(x,t)  is 
separable  in  space  and  time.  Therefore  a  possible  solution 
of  this  boundary  value  problem  is  of  the  form 


y(x,t)  =  Y  (x)  P  (t) 


(10) 


where  Y(x)  represents  the  general  beam  configuration  in 
vertical  vibration  and  depends  on  x  alone,  and  where  F(t) 
indicates  the  type  of  motion  the  beam  executes  and  depends 
on  t  alone.  Introducing  (10)  into  (9) 


-1 

Y  (x)  m(x) 


fEKxJ—^^FTtT 

3x 


a2F(t) 

3t^ 


(11) 
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where  partial  deriviatives  have  been  replaced  by  total 
derivatives.  In  addition,  the  variables  have  been  separated 
so  that  the  left  hand  side  depends  on  x  along  and  the  right 
hand  side  on  t  alone.  It  is  apparent,  therefore,  that  both 
sides  must  equal  a  constant.  Let  X  be  that  constant,  then 

d-  ^ i t-  -  XF (t)  *  0  (12) 

dtz 

where,  since  Y(x),  m(x)  and  I(x)  are  real,  X  is  real. 

Let  F ( t)  be  of  exponential  form 

F  ( t)  =»  Aest  (13) 

then  substitution  of  (13)  into  (12)  yields 
As2e3t  -  XAest  *  0 

S  t 

or  on  dividing  by  Ae  (for  Ae  /  0) 
s2  -  X  *  0 


so 


s  *  ±/  X 


(14) 


i ! 
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If  X  is  a  positive  number  then  the  two  roots  (14)  are  both 

real  and  equal  but  opposite  in  sign.  This  suggests  that 

there  are  two  solutions  for  F(t),  one  increasing  and  one 

decreasing  exponentially  with  time.  Both  of  these  solutions 

are  inconsistent  with  an  undamped  conservative  system  so 

the  possibility  of  X  being  positive  can  be  discarded. 

2 

Letting  X  =  -gj  we  have  s  =  ±ia  and 


F(t) 


„  iat  ^  .  -iat 
A^e  +  A2e 


(15) 


Expansion  of  this  expression  where 

i  tut  . 

e  =■  cos  cut  +  x  sxn  at 


yields 


F(t)  =  A^cos  at  +  iA^sin  at  +  A2COS  at  -  iA2Sin  at 
*  (A^+Aj )  cos  at  +  KA^-AjJsin  at 
This  must  reduce  to  a  real  expression  and  will  if 


A^  +  A2  *  Ccos<{> 
ifA^-Aj)"  Csin$ 
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then 


F  ( t)  =  Ccos(ojt  -  <J>)  (16) 

Proceeding  with  the  other  part  of  (11)  yields 
2  2 

~~y  [El  (x) d  Y4— )  =  w2  m(x)  Y(x) 
dx  dx^ 

0  <  x  <  L 

No  general  closed  form  solution  of  (17)  exists.  How¬ 
ever,  for  certain  special  cases  solutions  can  be  obtained  and 
for  illustrative  purposes  the  following  major  simplification 
will  be  made,  EI(x)  and  M(x)  are  made  constant  along  the 
beam's  length.  This  takes  all  ships  with  the  possible  ex¬ 
ception  of  some  box  shaped  barges  out  of  consideration  but 
the  information  derived  from  such  a  step  will  provide  some 
weight  into  the  problem.  (17)  reduces  under  these  assumptions 
to 


i 


(17) 


d-4J.^*I  -  34Y(x)  =0  S4  =  Ij2  (18) 

dx4 

The  general  solution  of  (18)  can  readily  be  verified  to  be 

Y (x) -A1sinsx+A2cosgx+A3sinhsx+A4coshsx  (19) 
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A2  =  A4  and  A1  *  A3  (22) 

At  x^L  substitution  into  (20)  yields 

A^ (sinhBL-singL) +A2 (coshSL-cosBL) *0  (23) 

and  substitution  into  (21)  yields 

A^ (coshSL-cosSL) +A2 (sinhgL+sinBL) =0  (24) 

Parenthetically,  these  two  equations  in  two  unknowns  can  be 
reduced  to  one  equation  leaving  A2  in  terms  of  A1 


19 


(25) 


-A^ (coshBL-cosBL) 
A2  =  sinhSL+sinBL 


Inserting  (25)  and  (22)  into  (19)  yields 

Y(x)»A[sin8xtsinhBx-  (g-j^|^“g°^)  (cosBx+coshBx)  ]  (26) 


From  linear  algebra  the  determinant  of  the  coefficients  in 
(23)  and  (24)  must  equal  zero  for  non- trivial  values  of  A^ 
and  A2  to  exist.  That  is 


sinhBL-sinBL 

coshBL-cosBL 


coshBL-cosBL 

sinhBL+sinSL 


0 


This  is  known  as  the  "characteristic  determinant"  and  reduces 
to 

2 

(sinhBL-sinBL) (sinhBL+sinBL) - (coshBL-cosBL)  *  0 


After  simplification  this  yields 

coshBLcosBL  =  1  (27) 

This  is  known  as  the  "characteristic  equation"  and  it  is 
apparent  that  it  can  be  solved  for  an  infinite  number  of 
values  of  3,  that  is, for  Sr  where  r*l,2,3...  which  are  defined 
by 


20 


are  obtained. 


From  (28)  an  infinite  set  of  frequencies,  ur. 

These  are  natural  frequencies.  For  each  ojr  equation  (26) 
yields  a  prescribed  vibration  shape  Yr(x)  known  as  a  "natural 
mode . " 

The  general  problem  of  determining  the  values  of  B  in 
(18)  for  which  nontrivial  solutions  Y(x)  exist  where  certain 
boundary  conditions  are  prescribed  is  called  the  "characteris¬ 
tic-value"  or  "eigen-value"  problem  with  the  parameters  ur 
being  the  eigenvalues  and  the  functions  Y^ (x)  being  eigen¬ 
functions.  [19] [25]  It  is  to  be  noted  that  the  eigenfunctions 
are  precise  as  to  configuration  but  not  to  absolute  value. 
Since  (18)  is  homogeneous  and  Yr(x)  is  a  solution,  then  so 
is  aYr(x)  where  a  is  a  constant  multiplier.  A  unique  property 
of  these  eigenfunctions  forms  the  basis  for  this  entire  train 
of  analytical  thought.  For  the  general  case  of  beam  vibration 
as  expressed  in  (6)  this  property  known  as  "orthogonality" 
is  expressed  by 

L 

/  m(x)  Yr(x)  Ys  (x)  =  0  r^s  (29) 

o 

This  orthogonolity  is  with  respect  to  the  mass  function  which 
serves  as  a  weighting  term. 


This  orthogonality  ii»p'.ies  that  any  possible  configura¬ 
tion  of  the  beam  can  be  represented  as  a  linear  combination 
of  the  natural  modes  of  the  beam.  That  is 


y<x)  =  2  C  Yfx) 
r»l  r  r 


This  is  true  for  any  and  all  time  t  so 


y(x,t)  *  Z  ri  (t)  Y  (x) 
r-1  r  r 


where  the  Cr  have  been  replaced  by  the  time  functions  nr(t) 
which  are  known  as  the  "natural  coordinates." 

As  previously  stated  the  eigenfunctions  are  not  deter¬ 
mined  in  regards  to  absolute  value.  A  convenient  normaliza¬ 
tion  scheme  is  such  that  the  eigenfunctions  satisfy 


/  m(x)  T.  (x)  Y  (x)  *6, 
o  r  s  3 


r  f s*l / 2 / 3 . . 


where  <5„.„  is  the  dirac  delta  function 
rs 


0  r^s 


1  r*s 


In  the  uniform  beam  example  insertion  of  one  eigenfunction 
and  eigenvalue  into  (18)  yields 
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Multiplying  both  sides  by  Yg (x)  and  integrating  over  the 
length  gives 


L  d4Y  (x)  L  2 

/  Y  (x)  El  - 5 -  =  /  oT  m  Ys  (x)  Yr  (x) 

o  s  dx*  o 


(34) 


Now  substituting  the  assumed  solution  (31)  into  the  beam 
equation  (9)  and  assuming  a  uniform  beam  yields 

at>  oo  d4Y  (x) 

2  ft  (t)  m  Y  (x)  +  I  q^tt)  El - 1 -  »  0  (35) 

r»l  r  r»l  r  dx* 

Multiplying  through  by  Yg (x)  and  integrating  over  the  length 
and  considering  (32)  and  (34)  yields  the  independent  set  of 
ordinary  differential  equations 

nr(t)  +  nr(t)  -  0  r-l,2,3...  (36) 

The  solution  to  each  of  these  equations  can  be  shown  to  be 


nr  *  Cr  cos(wrt  -  $r)  (37) 

where  Cr  and  $r  are  determined  from  initial  conditions  of  the 
beam,  y(x,o)  and  y(x,o). 
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Using  (37)  and  (j!)  and  setting  t=0  yields 


y(x,o)  =  2  C  cos<?  Y  (x) 
r=l  r  r  r 


y(x,o)  »  I  tC  u  sin*  Y  (x) 
r=l  r  r  r  r 


(38) 


multiplying  by  mYs(x)  and  integrating  over  the  length  yields 

L 


C  cosO  =  /  mY_.(x)  y(x,o)  dx 
r  0 

L  _ 

CrsinOr  =  /  —  Yr(x)  y (x,o)  dx 


(39) 


which  can  be  solved  for  Cr  and  4>r» 

The  free  vibration  case  for  the  uniform  free-free 
beam  has  now  been  completely  solved.  The  natural  modes 
Yr(x)  and  the  natural  frequencies  u)r  have  been  computed. 

It  is  the  orthogonality  of  these  modes  and  the  fact  that 
they  span  the  vector  space  in  which  they  are  contained  that 
makes  the  expansion  theorems  (30)  and  (31)  correct.  The 
natural  coordinates  nr(t)  can  be  determined  from  initial 
conditions . 

The  forced  vibration  case  is  slightly  more  complex. 
Here  substitution  of  the  assumed  solution  (31)  into  the 
beam  equation  (6)  yields 
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00  00  d4Y  (x) 

Z  H  (t)raYr(x)+Z  q  ( t)  El  - =-7 -  *  f(x,t)  (40) 

r=l  r  r=l  r  dx 

By  multiplying  by  Y^fx)  and  integrating  over  the  length  and 
considering  the  normalization  integrals  (32)  and  (34)  the 
following  independent  set  of  ordinary  differential  equations 
are  obtained 

nr(t)  +  uijnr(t)  =  Fr(t)  (41) 


where 


L 

F  (t)  *  /  Y  (x)  f  (x  / 1)  dx  (42) 

0 

The  solution  to  (41)  generally  involves  the  use  of  the 
convolution  integral 

1  ** 

n„(t)  =  —  /  F  (t)  sin  u  (t-r)dx 

r  wr  0  r  r 

+  n  (0)  cos  w  t  +  f|  (0)  sin  aj  t  (43) 

r  r  r  r 

where  n  (0)  and  n  (0)  are  as  determined  by  (39). 
r  r 

This  completes  the  classical  analysis  of  the  forced 
vibration  of  the  free-free  uniform  beam.  It  is  not  difficult 
to  see  that  significant  complications  arise  when  the  desire 
is  present  for  a  closer  representation  of  reality.  In  (7) 
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A  free  vibration  analysis  would  proceed  with  the  solution  of 
(9)  as  before  but  no  assumption  as  to  uniformity  would  be 
made.  The  resulting  general  solution 

<Mx)  *  AF(x)  +  BG (x)  +  CA(x)  +  DJ(x)  (44) 


can  be  obtained  for  sufficiently  well  behaved  SI(x)  and  m(x). 
Boundary  conditions  (7)  and  (8)  require  that 
AF" (0)  +  BG" (0)  +  CH"  (0)  +  DJ" (0)  =  0 

AF" (L)  +  BG” (L)  +  CH" (L)  +  DJ" (L)  =  0 

AF*(0)  +  BG,"(0)  +  CH'"(0)  +  DJ"'(  0)  =  0 

AF^L)  +  BG"*{L)  +  CH"'(L)  +  DJ"'(L)  »  0 
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This  leads  to  the  characteristic  determinant. 


F"  (0) 

G"  (0) 

H"  (0) 

J*‘  (0) 

F"  (L) 

GH(L) 

H”  (L) 

JM(L) 

F"*(  0) 

G'"(0) 

H'"(0) 

J’"(0) 

F"*(L) 

G'"(L) 

Hm(L) 

JM,(L) 

This  is  an  equation  in  w  and  its  roots  are  the  system 
natural  frequencies.  For  each  root  Lor  there  correspond 
specific  values  of  A,  B,  C,  and  D  and  therefore  a  specific 
eigenfunction  to.  The  nature  of  the  four  solutions  F,  G,  H, 
and  J  may  make  the  solution  of  the  determinant  quite  complex. 

This  is  indeed  an  understatement.  The  general  solution 
(44)  will  be  a  combination  of  transcendental  functions  and 
infinite  series.  Solution  of  the  determinant  would  be  very 
time  consuming  and  costly.  Additionally,  inherent  in  the 
previous  work  has  been  the  assumption  that  shear  defection 
and  rotatory  inertia  are  small  factors.  This  is  reasonable 
for  the  lowest  vibration  modes  but  increasingly  becomes  a 
source  of  error  as  the  mode  number  increases,  say  for  three 
or  more  nodes.  The  shear  defection,  for  example,  is  about 
15%  of  the  bending  deflection  for  the  lowest  vibration  mode, 
that  of  two  nodes  for  L/D  *  15.  [14]  For  vibration  of  about 
four  nodes  shear  deflection  becomes  greater  than  bending 
deflection. 
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B.  Prohl  Sequence 

A  method  which  allows  account  to  be  taken  of  shear  and 
rotary  inertia  and  can  be  solved  for  the  eigenvalues  and 
functions  is  the  technique  devised  by  M.A.  Prohl  in  1945. 

It  was  originally  produced  for  calculating  the  angular 
velocities  at  which  flexible  rotors  would  resonate.  The 
actual  rotor  was  simulated  by  a  "lumped  mass"  system  in 
which  the  rotor  mass  was  concentrated  at  a  number  of  stations. 
Station  spacing  was  a  matter  of  convenience  though  obviously 
the  greater  the  number  of  stations  the  greater  the  accuracy. 
The  stations  were  connected  by  weightless  rods  whose  bending 
stiffness  and  shear  characteristics  were  identical  to  those 
of  the  rotor  section.  This  effective  discretization  allows 
problem  solution  by  numerical  techniques.  The  procedure  is 
directly  applicable  to  the  analysis  of  ship  hull  vibration 
as  will  be  shown.  Let 

nu  =  station  mass  plus  added  mass 

*  bending  moment 
Vi  =  shear  force 

=  vertical  deflection 
9  ^  =  slope 

*  length  between  adjacent  stations 
to  =  angular  velocity 

KGA^  *  shear  stiffness  which  is  constant  along  L^ 

EI^  *  bending  stiffness  which  is  constant  along 


■  mass  moment  of  inertia 

and  i  3  subscript  denoting  consecutive  stations. 

The  equations  from  strength  of  materials  which  govern  the 


fV 


motion  of  the  stations  and  which  include  shear  and  rotary 
inertia  terms  are 


M.L. 


_  „  .  *i"i+l  ,  Vi+lLi+l2 

9  i+1  9i  +  El.  . 7  "2eT - 


'i+1 


'i+1 


Mi+1  *  Mi  +  vi+lLi+l  -  Ji+1  “  9  i+1 


v  =  V  +  9  L  ,  ^i+l2  ,  Wl+l3 
yi+l  yi  9iLi+l  2EI, 6EI, 


(1) 


(2) 


Vi+lLi+l 


i+1 


’i+1 


KGA 


i+1 


V 


i+1  *  Vi  +  mi“  yi 


(3) 

(4) 


It  is  to  be  noted  that  the  mass  term  nu  includes  the  mass  of 
hull  and  cargo  and  the  entrained  water  or  "added  mass." 

This  will  be  discussed  in  Section  C. 

The  relationships  between  the  loading  (shear) ,  bending 
moment  and  deflection  curves  are  well  known  from  strength  of 
materials.  As  u>  increases,  the  effect  of  including  the  rotary 
inertia  and  shear  deflection  terms  in  the  equations  becomes 
obviousr  the  shear  and  moment  terms  are  altered  increasingly. 

The  boundary  conditions  for  this  free-free  beam  are 
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VS  -  VB 


where  subscripts  s  and  B  represent  the  stern  and  bow  respec¬ 
tively  (more  accurately  the  after  terminus  and  the  forward 
terminus) .  Thus  for  station  one 


MSL1  .  vA 

ei1  IeT^ 


M  +  V1L1-J,oj28.  »  -Joj2  9 
s  1111  s 

ML2  V_L?  V. 

ys  +  9sLl  +  2EI^  +  GEI^  KGA^ 


yl  *  ys  +  9sL] 


It  can  be  seen  by  inspection  of  equation  (1)  through 

(4)  that  the  values  for  the  bending  moment  and  shear  force  at 

any  point  in  the  ship  will  be  a  function  of  y  ,  9  and  oi  in 

s  s 

addition  to  the  calculated  ship  parameters.  If,  at  the  start, 
y  is  set  equal  to  one,  a  unit  defection  of  the  stern,  then 

9 

the  equations  for  the  shear  and  bending  moment  at  the  bow  will 


be  expressed  by  the  following: 


j 


VB  -  a  +  b8s 
«B  *  -e  +  <S9s 

where  a,  b,  c  and  d  may  be  functions  of  w.  From  the  boundary 

condition  requirement  it  is  known  that  both  of  these  equations 

be 

must  equal  zero.  Arbitrarily  setting  Mg  *  0  and  Vg  ■  a  + 

and  performing  the  sequence  calculation  a  number  of  times  for 

various  values  of  u  a  curve  of  V.  versus  u  can  be  obtained. 

0 

This  curve  represents  the  exciting  force  V_  which  would  need 

O 

to  be  applied  at  a  specific  frequency  u>  to  generate  a  unit 
deflection  at  the  stern.  The  zeros  of  this  curve,  that  is 
where  the  curve  crosses  the  w-axis,  represents  the  natural 
frequencies  of  the  vessel  where,  with  the  absense  of  damping, 
a  vanishingly  small  exciting  force  would  set  up  synchronous 
motion.  With  the  natural  frequencies  obtained  the  governing 
equations  (1)  through  (4)  will  yield  the  displacement,  shear, 
moment  and  slope  for  each  station.  The  displacements  for 
each  frequency  uir  are  considered  for  this  discrete  case  to 
be  equivalent  to  the  continuous  system  eigenfunctions  and 
are  called  eigenvectors.  Indeed  as  the  number  of  stations 
becomes  large  the  eigenvector  will  converge  to  the  eigenfunc¬ 
tion.  As  for  the  continuous  case,  the  eigenvector  must  be 
normalized.  A  convenient  normalization  scheme  is  according 
to  the  following 
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1 


(9) 


N 

Z 

i=l 


which  is  the  discrete  representation  of  (32)  in  Section  A. 


C.  Added  Mass 

The  preceding  technical  treatment  of  the  free  and 
forced  vibration  of  a  ship  like  beam  involved  a  mass  term 
which  was  either  uniform  or  not  along  the  length.  This  term 
included  not  only  the  mass  of  ship  structure,  machinery, 
cargo,  etcetra  as  might  be  obtained  from  a  weight  curve  but 
also  a  hydrodynamic  "added  mass."  When  a  vibrating  body  is 
immersed  in  a  fluid  (air  or  water,  for  example)  the  fluid 
surrounding  the  body  is  put  into  continual  motion  which 
requires  the  expenditure  of  energy.  Under  conditions  of 
synchronous  or  resonant  motion,  the  body  with  this  entrained 
water  vibrates.  This  additional  mass  serves  to  lower  signi¬ 
ficantly  the  natural  frequencies  of  the  vessel.  It  is 
therefore  critical  that  a  determination  of  this  mass  be  made 
and  added  to  that  of  the  vessel.  The  total  mass  then,  called 
"virtual"  mass  is  the  sum  of  the  ship  mass  and  the  added  mass. 
This  phenomenon  has  been  observed  and  treated  for  a  few 
objects  in  classical  hydrodynamics.  In  the  case  of  an 
infinitely  long  circular  cylinder  having  a  cross-section  of 
radius  r  translating  perpendicularly  to  its  axis  with  velocity 
V  in  an  infinite  inviscid  fluid  it  is  possible  to  derive  the 


l 
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potential  flow  and  to  calculate  the  kinetic  energy  in  the 
fluid.  It  is  found  to  be  [20] 


T  *  ipr^V2 

per  unit  length  of  cylinder  where  p  is  the  fluid  density. 

2 

Letting  M'  *  pirr  it  is  apparent  that  M'  is  the  mass  of 
fluid  that  is  displaced  by  unit  length  of  the  cylinder.  It 
is  also  apparent  that  the  effect  on  the  cylinder  motion  from 
the  presence  of  the  fluid  can  be  allowed  for  by  an  addition 
of  M'  to  the  mass  per  unit  length  of  the  cylinder.  An  energy 
balance  equation  for  forced  motion  becomes 

^(|mv2  +  §M'V2)  =  FV 

where  F  is  the  external  force. 

On  simplifying 

dv 

(M+M')r~  =>  F 


or 


Mdv 

"dt 


F 
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Here  the  presence  of  the  fluid  is  equivalent  to  a  force  per 
unit  length  opposing  motion.  Obviously  this  virtual  ’'•ass 


effect  is  only  present  in  the  case  of  accelerated  mo cions 

.  v 

3uch  as  vibration.  For  a  sphere  of  radius  r  the  kinetic 
energy  of  the  fluid  has  been  shown  to  be  [20] 


T 


ih 


rV 


2  3 

Here  M'  =  jpirr  which  is  one  half  of  the  mass  of  displaced 
fluid.  In  general 


M'  =*  C  times  the  mass  of  displaced  fluid 


where  C  is  a  vertical  added  mass  coefficient.  Comparison 
of  the  added  mass  coefficients  for  the  cylinder  and  sphere 
both  of  radius  r  shows  the  former  to  be  twice  the  latter. 

This  reveals  the  three  dimensional  effect.  The  general 
method  for  determining  the  added  mass  coefficients  of  a 
ship  hull  is  a  stepped  procedure  utilizing  first  a  two 
dimensional  analysis  then  multiplying  by  a  three  dimensional 
correction  called  a  "J"  factor  which  will  be  discussed  later. 
The  discrete  analysis  of  the  Prohl  sequence  is  inherently 
compatible  with  this.  At  each  station  the  geometric  and 
material  characteristics  are  assumed  to  be  constant  along 
the  station  length.  Thus  the  station  can  be  considered  to  . 


be  a  portion  of  an  infinitely  long  cylinder  and  the  added 
mass  can  be  determined.  The  cylinder  is  supposed  to  consist 
of  the  underwater  hull  cross  section  and  its  above  water 
mirror  image.  The  actual  ship  added  mass  is  thus  one- half 
of  this  computed  amount.  This  analysis  is  based  on  the  fact 
that  the  fluid  motion  is  streamline  throughout,  that  is, 
irrotational .  This  is  reasonable  since  the  amplitudes  and 
absolute  velocities  are  small. 

This  assumption  of  irrotational  flow  implies  that  the 
curl  of  the  fluid  particle  velocity  vector  field  is  zero,  that 
is, 

V  x  V  *  0 

where  V  is  the  "del"  operator.  This  implies  that  the 
velocity  vector  field  is  the  gradient  of  a  potential 
function 

Y  a  V(j)  ^  ^  ) 

Additionally  where  the  flow  is  incompressible  and  there  are 
no  sources  or  sinks  a  fluid  continuity  equation  can  be 


written  in  the  form 


or  upon  substitution  of  (1)  into  (2) 

V2( j>  «  0  (3) 

This  is  the  well  known  Laplace's  Equation.  The  kinetic 
energy  of  any  flow  where  a  potential  exists  is  [20] 

2T  *  /  p(V<j>)2dV 

where  dV  is  a  volume  element  and  p  is  the  fluid  density.  As 
has  been  shown  the  added  mass  can  be  determined  from  the 
kinetic  energy  of  the  flow.  It  remains  then  to  compute  the 
potential  for  the  flow  around  a  ship  section. 

Classical  hydrodynamics  has  evaluated  the  potential 
for  a  very  few  geometric  shapes.  In  two  dimensions  flows 
past  circles  and  ellipses  have  closed  form  equations  and 
spheres  and  ellipsoids  of  revolution  have  equations  for  the 
potential  in  three  dimensions.  A  method  which  greatly 
expands  the  number  of  shapes  which  can  be  analyzed  is  called 
"conformal  mapping"  which  involves  application  of  analytic 
function  theory  and  functions  of  a  complex  variable. 


To  illustrate,  let 


where 


w  =  u  +  iv 

2  =  X  +  iy 

(4)  then  represents  a  mapping  of  all  points  in  the  z  plane 
onto  the  w  plane.  If  the  function  f(z)  is  single  valued 
then  the  mapping  is  one  to  one,  that  is  to  each  point  2 
where  f(z)  is  defined  there  exists  one  and  only  one  value 
of  w  =  u  +  iy  in  the  w  plane.  The  inverse  will  also  be 
true  if  f(z)  is  analytic  and  f '  (z)  ^  0.  If  these  conditions 
are  met,  then  it  is  also  true  that  relative  angle  and  shape 
are  preserved  in  such  a  mapping  (thus  the  term  name  "con¬ 
formal")  .  Right  angles  in  the  z  plane  will  be  mapped  to 
right  angles  in  the  w  plane.  Another  property  of  complex 
analytic  functions  is  that  the  real  and  imaginary  parts 
satisfy  Laplace's  equation  and  therefore  can  represent  fluid 
flows.  Knowing  that  the  potential  and  stream  functions 
satisfy  the  Cauchy- Rieraann  equations 


and  Laplace's  equation  (3)  an  analytic  complex  function  can 
be  generated  by  setting  the  potential  equal  to  the  real  part 
and  the  stream  function  equal  to  the  imaginary  part 

4>(2)  =  < p(x,y)  +•  ii p(x,y) 

This  function  is  often  called  the  complex  potential.  Now 
returning  to  (4)  and  writing  its  inverse 

2  =  F  (w) 

it  is  apparent  that  the  complex  potential  $>(z)  is  expressible 
as  a  function  of  w  in  the  form  $ [f (w) ] 

Further  this  function  is  analytic  and  therefore  repre¬ 
sents  a  new  fluid  flow  in  the  w  plane.  Each  of  the  equipo- 
tential  lines,  boundary  lines  and  streamlines  in  the  original 
plane  are  mapped  into  the  w  plane.  The  orthogonal  relationship 
is  present  and  thus  the  actual  flow  around  the  transformed 
boundary  is  known.  For  example,  the  flow  past  a  two  dimensional 
body  of  arbitary  shape  can  be  determined  (theoretically)  by 
generating  a  transformation  function  w  -  f(z )  which  will  map 
the  points  of  the  boundary  of  a  unit  circle  in  the  z  plane 
onto  the  boundary  of  the  arbitary  shape  in  the  w  plane.  The 
inverse  of  the  function  is  then  substituted  into  the  complex 


potential  for  the  flow  in  the  2  plane  which  is  well  known. 

Thus  a  complete  solution  for  the  flow  about  the  arbitary 
shape  is  obtained.  In  practice  the  inverse  may  not  be 
manageable  or  only  a  rough  approximation  to  the  desired 
shape  can  be  transformed  from  the  unit  circle.  Cross- 
sections  of  ships  can,  with  good  accuracy,  be  transformed 
from  a  unit  circle  however. 

The  first  and  definitive  work  in  this  area  was  pre¬ 
sented  in  1929  by  Professor  Prank  M.  Lewis  in  the  Transactions 
of  the  Society  of  Naval  Architects  and  Marine  Engineers.  The 
language  of  the  paper  is  a  bit  quaint  and  some  of  the  ques¬ 
tions  raised  have  long  since  been  answered,  but  the  method 
presented  for  computation  of  added  mass  has  seen  only  very 
minor  modifications  to  the  present  day. 

Let  x,y  be  the  coordinates  of  any  point  in  a  two 
dimensional  flow  and  (p,ip  the  values  for  the  velocity  potential 
and  stream  function  at  the  same  point.  The  equation 

2  =  x  +  iy  *  f  ($  +  i’i» ) 

where  f  designates  an  arbitrary  function  represents  a  two 
dimensional  irrotational  flow.  The  flow  in  the  x,  y  plane 
can  be  further  transformed  to  any  other  plane  X,  Y  by  the 
relation 


Z  *  X  +  iY  *  f  1  (x  +  iy) 
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is  another  arbitratry  function.  Z  is  of  course 


where  f ^ 
still  a  function  of  <{>  +  ity. 

Take  the  flow  past  a  unit  circle 


or  in  polar  form 

r=l 


and  transform  it  by  means  of  the  relation 
Z  *  X+iY  ■  Z+f+%  *  z+az“1+bz"3 


Now 


z 


iQ 

x  +  iy  =  re 


(5) 


(6) 


so 


x  +  iY 


ei9  +  ae-i9  +  be“ei0 


(7) 


Since 


*  cos8  +  isin9 
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it  follows  on  equating  real  and  imaginary  parts  that 


X  =  (1+a)  cos0  +  bcos38 

Y  =  (l-a)  sin6  -  bsin39 


(8) 


These  are  the  parametric  equations  of  the  transformed  circle. 
The  semi-axis  parallel  to  the  flow  (6=0)  will  be  of  length 


=  1  +  a  +  b 

and  the  semi-axis  perpendicular  to  flow  (9=ir/2)  is 


Y^  =  1  -  a  +  b 

so  the  ratio  of  half  beam  over  draft  will  be 


H  =  B/D 


1-a+b 

1+a+b 


(9) 


For  the  flow  past  a  unit  circular  cylinder  p  and  'p  are  given 
by 

$  *  U(r  +  i)  cosQ  ij/  =  U(r  -  p)  sinQ 


and  on  the  boundary  n»l 


<J>  =  2UCOS9 


V  *  0 


(10) 


This  represents  the  flow  at  the  boundary  of  a  stationary 
cylinder  and,  also  represents  the  flow  at  the  boundary  of 
the  transformed  cylinder.  To  obtain  the  solution  for  a 
moving  cylinder  in  a  fluid  stationary  at  infinity  a  super¬ 
position  of  a  uniform  flow  of  -U  must  be  made.  For  this 
uniform  flow 


<p  =  -UX  and  \p  *  -UY 


so  using  (11)  and  (8)  in  (10) 


<{>  *  U[(l-a)cos8  -  bcos38] 


ip  *  -U[(l-a)sin9  -  bsin39] 


The  kinetic  energy  of  the  fluid  is 


2T  *  p//(V2  +  V2)dxdy 

a  y 


and  since 


v  .  2i  v  *  — 

x  3x  y  3y 


2T  *  +  (|^)2]^y 


(ii) 


(12) 

(13) 
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An  application  of  Green's  theorem  allows  the  double 
integral  to  be  replaced  by  a  line  integral  computed  on  a 
path  which  encloses  the  area.  In  this  case  the  fluid  is 
bounded  by  the  cylinder  surface  and  by  a  circle  whose 
radius  is  allowed  to  approach  infinity. 


2T 


p/yn||)2  ♦ 


<||)2]dxdy 


-  P/  *|A  ds 
Cl+C2 


(14) 


where  ds  is  an  element  of  path  length,  3n  is  an  inwardly 
pointing  element  normal  to  the  curve,  is  the  cylinder 
circumference  and  C 2  represents  the  circle  with  radius 
approaching  infinity.  Since  the  fluid  is  at  rest  far  from 
the  cylinder  the  line  integral  over  C2  will  be  zero.  Since 

3  <b  3d) 

3n  *  ”  Ts 

the  kinetic  energy  is  given  by 

2T  =  pf<pdtp  (15) 

C1 

which  considering  that  the  expressions  for  $  and  <p  are 
dependent  only  on  3  equals 

8»2tt 

2T  ■  of  U [ (1-a) cos9-bcos38] [-U[ (1-a) cos0-3bcos38] ]d3 
0 

-  ttU2p  [(l-a)2+3b2]  (16) 
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If  the  semi-axis  of  the  transformed  circle  is  to  be  B 
instead  of  (1-a+b)  this  energy  becomes 


2T  -  iru2p[(l-a)2+3b2] 


B 


(1-a+b) 


-  T02pB2t(l_-a)Wi 
(1-a+b) 


(17) 


As  previously  mentioned  the  kinetic  energy  of  flow  past  a 
circular  cylinder  of  radius  B  is 


2T 


2  2 
iru  pB 


M'U 


therefore  the  virtual  mass  coefficient  for  the  transformed 


circle  is 


J.A*  U 


For  the  ship  hull  floating  on  the  surface  this  is  reduced 
by  one-half  so 


M'  »  |cttB2p 


Solving  (18)  and  (9)  simultaneously  for  a  and  b  yields 


2(CH2-H2-H-1) 


(1-H) (1+b) 

(I+h) 


and  the  parametric  equation  of  the  transformed  circles  when 
of  draft  D  instead  of  (1+a+b)  are 


„ r ( 1+a) cos8+bcos38 , 
Ul  (1+a+b)  J 


^  r  (1-a) sin0-bsin391 
Ul  (1+a+b)  J 


X  in  this  case  being  vertical  and  the  line  9=0  coinciding 
with  the  X  axis.  In  [22]  is  a  collection  of  graphs  for 
various  values  of  H(half  beam  to  draft  ratio)  and  C.  The 
H  values  correspond  to  common  ship  cross  sections.  For  each 
of  the  paired  values  of  C  and  H  the  constants  a  and  b  were 
computed  from  (21)  and  (20)  then  the  shapes  of  the  transformed 
circle  were  plotted  from  (22)  and  (23)  for  values  of  9  from 
0  to  t/2.  The  added  mass  for  any  ship  section  is  then 


1 


computed  as  follows:  a  visual  comparison  of  the  shapes 
in  [22]  is  made  with  the  body  plan  for  the  station.  Of 


course/  the  body  plan  must  be  properly  dimensioned  and  the 
value  of  H  known.  The  C  value  for  the  actual  ship  section 
is  obtained  by  interpolation  between  adjacent  curves  and  the 
added  mass  per  foot  for  the  ship  section  is  obtained  from 
(15)  remembering  that  B  is  half  beam.  This  procedure  is  the 
basis  for  present  day  calculations.  Some  refinements  which 
make  the  procedure  programmable  are  presented  in  a  paper  by 
L.  Landweber  and  M.  Macagno  in  1967  in  the  Journal  of  Ship 
Research.  An  equation  of  the  form 


(24) 


(of  which  (1)  is  a  truncated  version)  can  describe  either 
the  mapping  of  a  hull  shape  to  a  circle  or  a  circle  to  a 
hull  shape.  The  constants  A^  will  of  course  be  different. 
Furthermore,  inversion  formulas  exist  which  enable  the  co¬ 
efficient  of  one  transformation  to  be  determined  if  the 
inverse  coefficients  are  known.  Let  (24)  represent  the 
mapping  of  a  wetted  hull- shaped  contour  consisting  cf  a 
wetted  hull  and  its  mirror  image  across  the  waterline  from 
the  z  plane  into  the  w  plane.  S  is  the  bounded  area  of  the 

hull  contour  and  r  is  the  radius  of  the  circle  in  the  w 

o 

plane.  An  expression  for  the  added  mass  coefficient  C  is  [21] 


cv  -  rit^o+V-fi 

o 


(25) 


Let  the  inverse  of  (20)  be  of  the  form 


T 


(26) 


A  method  of  determining  the  coefficients  is  based  on  the 
property  that,  among  the  closed  curves  obtained  from  (22) 
by  varying  the  coefficients,  the  curve  which  will  bound  the 
maximum  area  will  be  a  circle.  Let  the  ship  curve  be 
termed  CT  and  its  mapping  in  the  t  plane  be  termed  If 

2  *  r(0)el9  then 

6  ~z  dz  *  6  r(dr+ird6)  *  i  6  r2de  =*  ZiS 
CT  CT  CT 

where  z  is  the  complex  conjagate  of  z.  Simlarly  for  C. 


dt  *  ZiS. 

C1 

where  S1  is  the  area  bounded  by  C^.  Also 


T 


(27) 


(28) 


and 


dT 


(1  - 


Bi 

7 


3B. 


.  .  )dz 


(29) 
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After  substitution  of  (28)  and  (29)  into  (27)  the  condition 
that  be  a  maximum  that  is 


x  —  1,2, .. .n 


yields  a  set  of  n  linear  equations  for  the  coefficients  , 
i  s  i  to  n.  Standard  matrix  manipulation  techniques  are 
available  for  then  determining  the  B^.  The  difficulties 
with  this  method  are  that  the  convergence  rate  with  respect 
to  n  is  relatively  low  and  the  accuracy  of  some  intermediate 
terms  decays  with  increasing  n.  The  solution  to  this  problem 
is  to  cut  off  the  first  mapping  at  n=8  and  then  map  the 
nearly  circular  curve  from  the  t  plane  into  a  much  more 
nearly  circular  curve  C2  in  the  w  plane  where 

C1  C2  C3 

w-t  +  ri  +  -4+-f.  .  . 

t  t 


The  are  determined  in  much  the  same  way  as  were  the  B^ . 

The  coefficients  of  equation  (24)  can  then  be  determined 
from  the  intermediate  coefficients  B^  and  C^.  The  added  mass 
coefficient  Cv  is  then  determined  from  (25) . 

This  thesis  makes  use  of  the  above  outlined  procedure 
to  compute  the  added  mass  for  vertical  vibration  of  hull 
sections  as  defined  by  offset  data  for  each  section. 
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ftp 


Another  significant  item  in  this  treatment  is  to 
calculate  the  three  dimensional  effect  mentioned  briefly 
earlier.  Lewis,  in  his  paper,  talks  of  a  "J"  factor  where 

T  _  Actual  Kinetic  Energy  of  Fluid 
w  “  Kinetic  energy  o£  fluid  if  motion 
is  confined  to  transverse  planes 

J  obviously  depends  on  the  underwater  configuration  and  will 
approach  unity  as  the  ratio  of  length  to  beam  gets  larger 
and  larger.  That  is  to  say  as  variations  of  the  hull  from 
section  to  section  diminishes. 

An  ellipsoid  of  revolution  is  one  of  the  very  few 
shapes  for  which  an  exact  solution  has  been  obtained  for  the 
three  dimensional  flow.  It  is  this  shape  that  is  used  for 
the  computation  of  the  J  factor  for  a  vibrating  ship.  As 
the  vibration  mode  goes  up,  and  the  number  of  nodes  increases, 
it  is  apparent  that  the  movement  of  water  becomes  more 
longitudinal  and  less  transverse  so  that  values,  where  i 
represents  the  i^*1  mode,  will  decrease.  The  calculation  of 
J2  and  was  made  in  the  paper  by  Lewis.  The  treatment  is 
similar  to  that  presented  for  the  two  dimensional  case  in 
that  an  equation  like  (1)  for  the  kinetic  energy  of  the  fluid 
i  used 


2T  -  ds 


(31) 
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where  in  this  case  ds  is  an  element  of  area.  The  solution 
of  this  equation  can  be  obtained  due  to  a  closed  form  solu¬ 
tion  to  Laplace's  equation  having  been  found  for  this  flow. 
That  is,  the  potential  <p  for  vibration  of  two  and  three  nodes 
can  be  substituted  into  (31)  and  values  for  the  actual  three 
dimensional  kinetic  energy  can  be  calculated.  The  denominator 
in  the  expression  for  J  (30)  is  computed  by  summing  (integra¬ 
ting)  the  values  obtained  from  a  two  dimensional  analysis  of 
cross  sections  of  the  body.  Lewis  tabled  the  values  of  J2 
and  for  a  variety  of  ratios  of  L/B.  The  values  are 
shown  in  Table  1. 

The  analysis  performed  by  Lewis  required  that  the 
vibrating  cross  section  remain  in  the  same  plane,  that  is, 
the  motion  results  from  pure  shear.  This  does  not  reflect 
reality  and  consequently  there  has  been  much  discussion  on 
the  subject  beginning  with  a  written  comment  of  Lewis'  paper 
by  his  peer  J.  Lockwood  Taylor  [22] [30]  and  continuing  to 
the  present  day. 

Current  procedure  [32]  involves  the  determination  of 
J  by  accurate  measurement  of  the  natural  vibration  frequencies 
of  a  given  body  in  air  and  in  water.  For  Lewis'  circular 
ellipsoid  of  revolution  a  simple  analysis  is  possible  because 
the  distribution  of  added  mass  is  directly  proportional  to 
the  distribution  of  actual  mass  ( Cv  *  constant  along  the 
entire  length) .  The  mass  distribution  is 
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>  a  irmia 


(32) 


«b  =  «2pb 

.where  r  is  the  radius  of  the  circular  cross  section  and  a 
function  of  longitudinal  position,  p  is  density  and  subscript 
B  refers  to  the  body.  The  corresponding  added  mass  is 

M'  =»  CvJ7rr2pw  (33) 

where  Cv=l  and  subscript  w  refers  to  the  water  in  which  the 
body  is  immersed.  J  here  is  considered  to  be  a  constant 
value  applied  at  each  cross  section  so  that  the  actual  three 
dimensional  added  mass  is  obtained  from  the  two  dimensional 
analysis  performed.  From  basic  vibration  theory  [22] 


u)  =  Sk7m~ 

n 


where  K  represents  stiffness  and  M  mass.  Put  another  way 


f  0  M 
,  a.  2  w 

(C 


where  f  represents  frequency  and  subscript  a  represents  air. 
This  equation  for  the  circular  ellipsoid  considered  is 
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TABLE  1 


LEWIS 

L/B  J  2 

■  K  - - - - - 


TOWNS IN  LEWIS  TOWNS IN 

J3  J3 


oo  1.000 

30 

22.4  .954 

16 

15.82  .922 

12.95  .897 

12 

10.02  .859 

9 

8.48  .829 

7.11  .789 

7 

5.83  .741 

5.07  .704 

5 


.947 

.883 

.837 

.764 

.683* 


.593 


1.000 

.932 

.890 

.859 

.808 

.770 

.725 

.670 

.629 


.  924 

.842 

.793 

.696 

.591 
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f  -5  Pn/ffr2dL+Jp/-n-r2dL 

.  a.  £  o  W 

fw  "  PB/7rr2dL 


»  1 


(35) 


where  the  vibrating  mass  in  air  M  is  assumed  to  be  due  only 

cl 

to  the  mass  of  the  body  and  dL  is  differential  length. 
Equation  (35)  is  reduced  to 


J  = 


Ur) 

w 


(36) 


Values  from  this  equation  obtained  by  R.L.  Townsin  [32]  are 

listed  in  Table  1  and  a  comparison  with  the  Lewis  values 

shows  ?  significant  decrease  reflecting  the  fact  that  no 

restrictions  or  assumptions  were  made  concerning  the  cause 

of  the  motion.  For  ship- like  vibrating  bodies  the  procedure 

is  somewhat  more  complex.  Frequencies  in  water  and  air  ave 

both  measured  experimentally  and  calculated.  For  the  water 

calculation  several  assumed  values  of  J  are  used  to  obtain 

a  plot  of  Jn  versus  f.  A  reduction  constant  obtained  from 

comparison  of  the  two  air  frequencies  is  used  to  correct 

the  measured  water  frequency.  This  abscissa  is  used  to 

obtain  the  J  ordinate  from  the  J  versus  f  plot.  Townsin 
n  n 

proposes  and  this  thesis  employs  the  following  empirical 
formula  for  Jn  as  a  function  of  B/L 
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J  »  1.02  -  3(1. 2-1/N) B/L  n*2, 7  (37) 

n 

This  formula  agrees  closely  with  data  through  J7. 

The  remaining  consideration  which  is  presented  deals 
with  a  "sectional  longitudinal  inertia  correction”  labelled, 
by  the  author,  j ^  (the  lower  case  to  reflect  a  relation  with 
the  Lewis  J  factor) .  In  most  determinations  of  actual  added 
mass  per  section  an  equation  of  the  form  of  (19)  is  consis¬ 
tently  used.  That  is 

rai  *  I0!*3!3  (38) 


and 


— —  *  J 
Em! 

i 


N 

M'  *  J  Z  m.  (39) 

i-1 


so 


Mi  ’  2ciBipJ 


(40) 
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where  iik  represents  the  2-D  section  added  mass,  and 
are  the  section  added  mass  coefficient  and  half  beam  respect¬ 
ively,  M'  is  the  3-D  added  mass  for  the  entire  body  and 
is  the  supposed  3-D  section  added  mass.  The  departure  from 
reality  in  (40)  is  apparent.  J  of  the  entire  body  represents 
a  fixed  fraction  which  when  multiplied  by  the  added  mass  of 
a  2-D  flow  gives  the  actual  3-D  added  mass  for  the  entire  body. 
There  are  some  sections  of  a  vibrating  ship  where  the  fluid 
motion  is  in  fact  very  close  to  two  dimensional  and  the  added 
mass  should  be  very  close  to  that  calculated  from  (38) .  At 
other  sections,  particularly  near  bow  and  stern,  the  flaw 
departs  quite  significantly  from  a  pure  two  dimensional  flow 
and  the  added  mass  would  be  less  than  that  obtained  from  (40) . 
To  take  into  account  these  aspects  the  author  proposes  to 
use  ji. 

Let  the  added  mass  at  section  i  for  a  2-D  analysis  be 

mi  *  {38) 

and  let  the  actual  3-D  added  mass  be 


(41) 


then 


N  N 

M'  =  2  M?  =2  m . j .  (42) 

•  i  1  •  .  11 

1*1  1=1 


The  requirement  that  equations  (39)  and  (42)  yield  the  same 
value  for  M'  requires  that 


N  N 

j  Z  m!  **  2  m.  j .  (43) 

i-L  1  i*l  1  1 


Let 


h  -  1  -  dib 


(44) 


where  d.^  is  a  difference  factor  computed  for  each  section  by 

N 

n  £iml  ”i 

di  *  - 5 -  +  — - 

N 

i  n  I  m! 


[ mi+l~nii-l '  +  i=l 


2N 


(45) 


and  b  is  a  normalization  constant  to  be  determined  by  substi¬ 
tution  of  (44)  into  (43) .  In  (45)  N  is  the  number  of  stations 
and  n  is  the  number  of  nodes  in  the  vibrating  mode.  As  the 
number  of  nodes  increase  three  dimensional  flow  increases  at 
all  stations  so  An  example  will  serve  to  illustrate. 

Assume  that  an  analysis  has  determined  the  2-0  sectional 
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added  masses  to  be  as  listed  in  column  2  of  Table  2.  The 
difference  factors  from  (45)  are  listed  in  column  3.  Sub¬ 
stituting  the  totals  of  columns  2  and  4  into  (43)  yield 

.81(692)  *  692  -  (28380  +  47886n)b 

b  »  131.5/(28380  +  47886n) 


b(xl0  3)  n 

4.633  0 
1.723  1 
1.059  2 
0.764  3 


Column  5  lists  or  the  m|J.  Columns  6,  7  and  8  show  the 
actual  added  mass  at  each  station  for  vibrations  of  0,  1  and 
3  nodes  respectively.  Notice  as  n  increases  the  added  mass 
approaches  the  values  in  column  5. 
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D.  Wave  Excitation 


The  forcing  function  f(x^,t)  used  in  the  previous 
treatment  represents  the  total  external  force  per  unit  length 
at  any  section  x^  and  any  time  t.  It  includes  the  many  effects 
of  the  waves  encountered,  the  motion  of  the  ship  and  the  fluid 
flow.  It  is,  as  such,  very  complicated  and  the  subject  of 
much  research.  [12] [16] 

The  most  simple  formulation  of  wave  excitation  would 
be  due  strictly  to  buoyancy  changes  from  wave  passage. 

V 

f  ( x . ,  t )  -  /  y  dV 
Vo 

where  an  elemental  volume  dV  is  equal  to  b^y  to  give 

C(t) 

f (x. ,t)  =  /  yb.dy  (1) 

1  0  i 

where  b^  is  the  breadth  of  the  ship  and  a  function  of  y, 

?  is  the  wave  height  at  any  time  t  and  y  is  the  specific 
weight  of  water.  It  has  been  traditional  in  naval  architec¬ 
ture  to  use  a  trochoidal  wave  [14]  to  model  the  wave  form 
for  deep  water  waves.  With  sharper  crests  and  flatter 
troughs  it  makes  a  very  close  approximation  to  the  actual 
sea  surface.  However,  this  wave  form  cannot  be  derived  from 
the  velocity  potential  for  surface  waves  and  the  parametric 
equations  which  describe  a  trochoid  are  an  additional 
complexity  in  any  analysis. 
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A  wave  form  which  is  more  manageable  and  is  derived 
from  hydrodynamic  theory  of  waves  of  finite  amplitude  was 
prepared  by  G.G.  Stokes  in  1847  [14]  . 

K  72 

C  *  SCOS(u)t-Kx)+^2—  COS2  (oit-Kx)  (2) 

where  ~Q  is  half  wave  height  from  trough  to  crest  and  the 
circular  frequency  oj  and  the  wave  number  K  are  related  by 

K  =  u2/g  *  2tt/Lw  »  g/V*  (3) 

for  deep  water  waves  where  L  is  the  wave  length  and  V  is 
the  celerity  or  wave  velocity.  In  (2)  a  simple  traveling 
cosine  wave  is  modified  by  a  harmonic  which  is  half  the 
length  of  the  fundamental  but  which  moves  with  the  same 
celerity.  This  wave  form  has  the  same  sharp  peaks  and  wider 
troughs  as  does  the  trochoid  but  the  fact  that  it  is  nonlinear 
with  respect  to  height  makes  it  somewhat  less  attractive.  In 
the  program  suitable  approximations  will  be  made. 

It  is  well  known  that  the  pressure  in  a  still  body  of 
water  is  equal  to  yh  where  h  is  the  distance  below  the  still 
surface.  It  is  also  well  known  that  in  deep  water  the  water 
particles  have  a  circular  motion  which  serves  to  reduce  the 
pressure  according  to  Bernoulli's  theorem.  That  is,  for  a 
simple  harmonic  wave. 


(4) 


Yh  + 


;ye“Khcos (wt-Kx) 


This  is  the  Smith  effect.  The  pressure  on  the  hull  due  to 
the  wave  passage  is  then 

p*  -^’Kh 


where  5  is  given  by  (2).  The  buoyancy  due  to  this  pressure 
is  an  integration  of  the  vertical  component  of  the  pressure 
over  the  surface  area  of  a  hull  slice  of  unit  length.  That 
is 


f  (x,t) 


s 

-/  PM  n-j  ds 
0  w 


where  n  is  an  outward  pointing  unit  normal  to  the  hull  and 
j  is  the  unit  vector  in  the  positive  y  direction.  For  a 
ship  vertically  sided  to  the  baseline  this  reduces  to 


f(x,t)  =  Pw(T)b 


or  with  (5) 


f  (x,t) 


b?ye 


-KT 


(6) 


where  T  is  the  draft. 
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Since  the  pressure  given  by  (5)  is  seen  to,  in  effect, 
reduce  the  waveheight  S,  an  effective  wave  height  s*  can  be 
defined  so  that 


fx(x,t)  *  ybC* 


(7) 


where  the  assumption  is  made  that  the  ship  is  vertically 
sided  at  the  waterline  so  that  b  is  constant  over  the  range 
of  s*.  The  subscript  1  indicates  that  this  buoyancy  force 
is  only  part  of  the  total  excitation.  is  given  by  [16] 


C*  * 


zeKydy] 


(8) 


where  z  is  the  hull  offset  and  y  is  distance  to  the  design 
waterline.  The  integral  serves  to  weigh  the  determination 
of  the  reduction  factor  by  the  hull  shape.  For  a  section 
with  a  vertical  side  to  the  baseline  (8)  would  reduce  to 

5*  -  ?[1-K  ?  eKydy] 

-T 

=  5(1 -/  eKyKdyJ 

-T 

=  5 [1- (l-e”KT) ] 
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so  that 


f^Cx/t)  =  YbSe 


which  is  identical  to  (6).  The  discrete  form  of  (8)  is 


S*  =  C[l-|  2  z.  eKyi  Ay  ] 
D  i=l  1  1 


and  for  the  discrete  forcing  function  per  unit  length 


(x. ,t)=  yb.?*=Yb.  ?[1-t—  2  z.  eKyi  Ay.] 
1  J  J  Dj  i»l  1  1 


where  ?  is  given  by  (2) . 

The  behavior  of  the  bracketed  term  is  to  reduce  the 
wave  height  as  K  increases.  For  small  K  the  term  is  very 
close  to  1.  Consider  the  straight  sided  vessel  again.  A 
reduction  in  wave  height  to  one-tenth  of  the  long  wave  value 
requires 


.  1  *  e 


-2.30  *  -KT 


K  »  T/2.30 


For  a  draft  of  20  feet  K  *  8.70  which  corresponds  to  a  wave 
length  of  .72  feet.  So  as  K  increases  the  correction  term 
approaches  zero.  Thus  there  are  two  major  effects  taking 
place  to  reduce  the  energy  available  for  excitation  of  the 
higher  modes  of  a  ship  hull.  First  the  wave  height  decreases 
as  the  frequency  increases  and  the  effect  on  an  immersed  body 
decreases  as  the  frequency  increases. 

Due  to  the  logarithmic  decrease  in  (9)  the  discrete 
form  of  the  correction  term,  hereinafter  called  D,  causes 
some  convergence  problems,  that  is, as  K  increases  D  fails  to 
approach  zero  if  the  increment  Ayi  is  too  large.  In  the  pro¬ 
gram  four  foot  waterlines  are  used  and  the  following  slightly 
modified  expressions  must  be  employed. 


D  - 


e 


-KTS 


(11) 


Equation  (11)  has  the  proper  limiting  values  and 
accounts  for  the  hull  shape  also. 

There  are  two  other  parts  to  the  forcing  function  to 
be  presented.  As  the  wave  passes  down  the  side  of  the  ship 
the  vertical  component  of  the  orbital  velocity  has  a  damping 
effect  and  this  force  can  be  expressed  by 


f _ (x • , t)  -  (N'-Vdm'/dx) ;* 
z  j 
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(12) 


m 


9 


where  (N'-Vdm'/dx)  is  the  sectional  hydrodynamic  damping 
coefficient  which  will  be  discussed  later  and  is  the 
time  derivative  of  the  effective  wave  height.  The  third 
force  is  due  to  the  inertia  of  the  added  mass 


f3(Xj,t)  *  IHj  ?* 


(13) 


where  %  is  the  effective  vertical  acceleration  of  the  wave 
height.  As  was  disctu  jed  in  Section  C  the  added  mass  may  be 
incorporated  into  the  equation  of  motion  in  two  ways,  either 
as  an  addition  to  the  mass  of  the  vessel  and  placed  on  the 
left  hand  side  or  as  an  exciting  term  and  placed  on  the  right 
hand  side.  The  superposition  principle  allows  the  separation 
of  the  added  mass  effect  in  this  case  since  two  separate 
motions  are  being  superposed,  the  vibratory  motion  of  the 
vessel  and  the  harmonic  motion  of  the  waves. 

In  (13)  nij  is  the  added  mass  calculated  for  immersion 
to  the  still  water  waterline.  As  the  wave  passes,  the  draft 
at  a  section  changes  and  so  does  the  added  mass,  but  here 
the  sides  at  the  waterline  are  assumed  to  be  nearly  vertical 
and  the  change  in  added  mass  is  very  small  so  the  effect  can 
be  neglected. 

The  total  excitation  force  is  then 


3 

f  (x,t)  *  2  f .  (x ,  t)  (14) 

i-1  1 
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and  the  modal  excitation  force  is  obtained  by  multiplying 
by  and  summing  over  the  ship  length  as  shown  in  Section  A. 

N 

F  (t)  =*  %  Y  (x.)f(x.,t)  (15) 

r  i=l  r  1  i 

The  correct  surface  height,  velocity  and  acceleration  are 

_2 

£*  =  D^cos  (ojt-Kx)  cos2(wt~Kx) 

£*  =  -Du^sin  (uit-Kx) -DKw£2  sin2(o>t-Kx)  (16) 

'£*  =  -Dw2  scos  (wt-Kx) -2DKa)2  cos2  (ait-Kx) 

and 

f(xi,t)  =Ybi;*+(N'-Vdm'/dx)  c*  +  m /** 

2— 

=  (ybi-m|Dia3  ;)  cos  (wt-Kx) 

(N '  -Vdm*  /dx)  sin  (ojt-Kx) 

Yb.D.K?2  2_2 

+  ( — i— i - 2D.Ku!  C  ni. )  cos2  (ut-Kx) 

M  1  • 

(N1  -Vdm'/dx)  sin2  (ut-Kx)  (17) 
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Let 


2- 


Gil  =  ^biDi^  '  miDiu  5 


Gi2  =  -DjW^N'-Vdm'/dx) 


T°iDil'.;2  2—2 

st3  =  ~ ¥ - ? 


Gi4  =  -D^uk2  (N*  -Vdm'/dx) 


then 


fCx^/t)  =*  G^coswtcosKx  +  siiiuitsinKx 


+  Gi2sinwtcosKx  -  G^coswtsinKx 


+  G^2cos2wtcos2Kx  +•  G^2sin2uTsin2Kx 


+  G^4sin2wtcos2Kx  -  G^4cos2uiTsin2Kx 


Substitution  into  (15)  yields 


N 

F  (t)  *  coswt  2  Y  (x.)  G . , cosKx . A 
r  i=1  r  i  il  i 

N 

+  sincot  l  Y  (x . )  G . ,  sinKx .  B 
i=l  r  1  11  L 
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where 


rl 


N  2  N  2 1 
=  {(Z  Y  (xi)GilcosKxi)  +  (Z  Yr  (xi)GilsiniCxi)  > 

i=l  r  i=l 


2  1 1/2 


N  2  N  2-.1/2 

5r2  =  {(Z  (xi)Gi2cosKxi)  +  (Z  (xi) Gi2sinKxi)  } 

“  (22) 

N  2  N  2-.1/2 

?r3  =  {(Z  Yr(xi)Gi3cos2Kxi)  +  (Z  ^  (x±)  G^sii^Kx^  } 


N  2  N  2 

5r4  =  {(S  iYr(xi)Gi4c°s2Kxi)  +(S  (xi)Gi4sin2Kxi)  } 


2,1/2 


J  Yr<xi}GilsinKxi 
9  ^  -*  arctan  N~ - - - 

I  Yr(*i)GilcosKxi 
1™* 


N 

Z  Yr(xi)Gi2cosKxi 
9  2  =  ir-arctan  - 


Z  Yr (x. )G. -sinKx. 

i_l  a  1.4  i 


N 


9 3  =  arctan  N 


Z  Yr(x. ) G. ,sin2Kx . 
i=l  1  _ _ 


Z  Y  (x . ) G . «cos2Kx . 
1=1  r  1 


(23) 


N 

Z  Y  (x.)G..cos2Kx. 

i=l  r  1  _ _ 

N 

?  ,Yr(Xi)Gi4  sin2Kxi 
i*l  r 


tt  -arctan 


The  Pr^  i-1,2,3,4  are  known  as  "participation  factors" 
and  were  obtained  by  weighting  the  wave  forces  by  the  mode 
shape.  The  factors  indicate  the  participation  of  the  wave 
forces  in  exciting  a  vertain  mode.  Solution  of  the  independent 
equations  for  the  modal  response  as  discussed  in  Section  C 
yields 


4 

n-(t)  *  z 

r  i=l 


Pric°s  (uit-8i+er) 

7T"2  2,2 “77  ,  2, 1/2 

(  (o)r-aj  )  +(Crw)  ) 


where 


£ 


r 


C  ui 

arc  tan  (— 5 — j-) 

01  -til* 


(24) 


The  most  significant  part  of  this  exact  time  history 
is  the  magnitude  of  the  response.  Equation  (24)  may  be  mani¬ 
pulated  to  yield  the  amplitude  of  the  sinusoidal  response 


n 


r 


{(2  P^cose.)2  +  (2  P^.sine. )  2}1/2 
i-1  ri  1  <-l  ri  1 


(25) 


The  nonlinearity  of  the  expression  for  the  height  of 
the  sea  surface  can  be  eliminated  by  dropping  the  modifying 
harmonic.  Although  the  simple  cosine  term  is  not  an  exact 
representation  of  the  sea  surface  it  is  ideal  for  representing 
the  many  component  waves  of  the  general  sea.  It  also  allows 
the  principle  of  superposition  to  be  applied  and  for  these 
reasons  the  program  encorporates  this  modification.  Result- 
antly  G^G^P 3»P 4=*0  with  significant  subsequent  simplification 
of  the  program. 
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This  analysis  assumes  that  the  presence  of  the  vessel 


does  not  disturb  the  wave  motion  which  it,  in  fact,  does. 

This,  then,  provides  an  upper  bound  to  excitation  and  real 
forces  will  be  less.  How  much  less  is  a  matter  for  speculation 
and  research. 

E.  Damping 

It  would  not  be  incorrect  to  say  that  relatively  little 
is  known  about  damping  associated  with  ship  hull  vibration. 

To  be  able  to  calculate  the  amplitudes  of  desired  quantities 
such  as  stress  levels,  bending  moments  and  displacements, 
the  distribution  of  damping  in  the  hull,  cargo  and  surrounding 
water  should  be  known.  The  available  literature  is  sparse, 

[4] [11] [16] [17]  but  it  is  generally  assumed  that  energy  is 
dissipated  thru  a  number  of  mechanisms.  A  convenient  cate¬ 
gorization  divides  the  mechanisms  into  those  which  are  hydro- 
dynamic  and  those  which  are  not.  Water  friction  and  the 
generation  of  both  surface  and  pressure  (sound)  waves  are 
the  hydrodynamic  effects.  Structural  and  cargo  damping  are 
non-hydrodynamic.  Another  possible  way  of  differentiating 
the  loss  mechanisms  is  to  determine  where  the  damping  takes 
place,  either  within  or  external  to  the  hull. 

As  described  previously,  superposition  of  the  natural 
modes  of  the  ship  yields  for  vertical  vibration 

00 

y (x, t)  »  Z  Y  (x)n  (t)  (1) 

r-0  r  r 


v-VVW 
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where  Y  (x)  is  the  normalized  mode  shape  and  n  (t)  is  the 

r  r 

natural  coordinate  which  incorporates  the  natural  frequency, 
initial  conditions  and  modal  damping  factor  into  a  time  de¬ 
pendent  function.  Motion  for  r=0,l  is  heave  and  pitch 
respectively.  Fortunately  the  damping  in  ship  hull  vibration 
is  relatively  small  so  that  calculations  for  the  natural 
modes  and  frequencies  can  be  carried  out  assuming  no  damping 
without  significant  error.  This  is  especially  true  for  the 
lower  modes. 

Other  qualitative  information  concerning  damping  has 
been  known  for  some  time.  Higher  modes  are  less  affected  by 
hydrodynamic  damping  due  to  in  part  the  smaller  displacements 
involved.  In  these  modes  structural  damping  is  the  major 
factor.  The  situation  is  reversed  for  lower  modes.  The 
generation  of  sound  waves  increases  as  the  vibration  fre¬ 
quency  goes  up  but  this  remains  a  very  small  quantity  through¬ 
out  the  range  of  interest,  that  is  oj  less  than  about  35  rad/ 
sec.  For  rigid  body  motion  there  is,  of  course,  no  structural 
damping.  If  vibration  analysis  is  of  concern  and  sea-keeping 
is  not  then  discussion  can  be  limited  to  distortion  modes 
r  -  2.  Surface  wave  generation  increases  with  decreasing  a>. 
For  even  the  largest  of  vessels,  however,  ^  is  greater  than 
about  2.5  rad/sec  and  wave  making  damping  is  far  less  than 
structural  damping. [17] 
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In  a  discussion  of  damping  it  is  common  to  refer  to 
the  logarithmic  decrement  <5  associated  with  a  mode  or  given 
damping  mechanism. 

t 

5  *  27T5//1-52  =  ln(x1/x2)  (2) 

where  c  is  the  critical  damping  ratio  associated  with  the  mode 
or  mechanism  and  x^  and  x ^  are  values  of  the  response  separated 
by  one  wave  period. 

Since  structural  damping  is  the  significant  factor  and 
appears  to  be  manageable  there  has  been  some  effort  directed 
at  its  analysis. [4]  This  damping  derives  from  hysteresis  in 
the  steelwork.  The  logarithmic  decrement  for  this  is  a  power 
function  of  the  stress.  For  levels  of  stress  less  than  about 
30,000  psi 

5  =  Co11  (3) 

where  C  is  a  material  constant  and  n  =  3 .  At  higher  stress 
the  exponent  increases  rapidly  and  other  effects  such  as 
stress  history  become  important.  For  steel  plate  (3)  becomes 

6  x  10“5  -  5  -  1.6  x  10-4 

i 
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Values  for  a  welded  structure  go  up  with  the  complexity  of 
the  structure  and  may  lie  in  the  range 

9  x  10"3  -  5  -  3  x  10"2 

This  dramatic  increase  is  explained  in  part  by  the  high 
residual  stresses  from  welding  and  stress  concentrations  at 
joints.  There  may  also  be  as  yet  undiscovered  hysteresis 
effects  which  are  structure  related  instead  of  material 
related. 

Ships  of  current  design  are  usually  welded  but  in 
riveted  ships  coulomb  friction  at  the  joints  is  very  important. 
In  traditional  analysis  this  dry  friction  is  treated  differ¬ 
ently  from  other  damping  effects.  The  coulomb  damping  force 
is  not  a  function  of  velocity  rather  it  has  a  constant  value 
for  all  velocities  and  varies  only  in  sign.  This  directional 
force  is  treated  as  an  excitation  force  which  opposes  motion 
and  equations  for  each  direction  arc  solved.  One  test 
concluded  that  riveted  ships  have  twice  the  structural 
damping  of  welded  ships. [17]  This  would  obviously  vary 
with  ship  age  and  history. 

Experiments  have  been  performed  on  loaded  and  unloaded 
vessels  to  determine  cargo  damping  effects.  Results  show 
that  cargo  does  increase  the  logarithmic  decrement  but  for¬ 
mulations  which  have  been  made  to  quantify  the  effect  have 
failed  to  predict  well. [17] 


The  literature  separates  the  generation  of  surface 
waves  from  the  other  damping  mechanisms.  It  is  that  part 
of  the  hydrodynamic  damping  that  can  be  derived  from  potential 
theory  and  the  only  part  of  ship  damping  with  a  solid  mathe¬ 
matical  base.  However,  the  literature  goes  on  to  say  that  in 
the  distortion  modes  the  natural  frequencies  are  relatively 
high  and  the  displacements  so  small  that  for  most  existing 
vessels  the  effect  can  be  taken  as  zero.  A  bit  ironic.  In 
equation  (11)  of  Section  D  the  term  in  parentheses  (N’-Vdirf/dx) 
is  the  sectional  hydrodynamic  damping  coefficient  according  to 
the  strip  theory  approach  used  in  this  paper.  N '  is  a  para¬ 
meter  which  quantifies  the  generation  of  surface  waves  and 
Vdirf/dx  is  a  forward  speed  correction  term.  V  is  velocity  and 
m'  is  the  added  mass.  This  coefficient,  it  is  seen,  varies 
along  the  hull  length  and  this  distribution  of  damping  effect 
is  important  for  computation  of  exciting  forces.  The  normal 
mode  analysis  in  Section  A  required  that  the  damping  be  pro¬ 
portional  to  the  mass  and  stiffness  distribution.  Otherwise 
the  modal  equations  would  not  be  uncoupled.  Since  the  pre¬ 
dominant  damping  source  is  structural  this  assumption  is  not 
too  far  off,  even  though  the  distribution  of  the  sectional 
hydrodynamic  damping  is  weighted  at  the  bow  and  stern.  N '  is 
nearly  zero  amidships  but  gets  larger  at  the  bow  and  stern. 

The  effect  is  still  very  small.  Only  the  forward  speed 
correction  will  be  considered  in  this  work  and  consequently 
the  amplitude  obtained  will  be  slightly  high. 


Consideration  of  equation  (2)  brings  up  the  fact  that 


in  much  practical  usage  all  the  damping  effects  are  lumped 
into  one  parameter  and  its  value  obtained  from  a  spectral 
analysis  of  stress  data  taken  from  a  vessel  which  is  loaded 
impulsively  so  that  all  modes  are  excited.  Several  full  scale 
measurements  from  slamming  or  mechanical  excitation  (anchor 
drop  test)  have  been  obtained  and  modal  damping  constants 
calculated. [11] [17]  Some  convenient  empirical  relations  have 
been  devised  from  such  experiments  and  other  means  so  that 
values  for  the  Cr  can  be  obtained  quickly.  From  basic  vibra¬ 
tion  theory  [24] 


2Vr 


- x - x- 

/ (2tt) 


(4) 


and  a  formula  for  the  logarithmic  decrement  which  seems  to 
give  intermediate  and  hopefully  representative  values  is 
[17] 


<$r  »  1.065  x  10"2  ur1/2  u)  <  31. 5  (5) 

Equations  (4)  and  (5)  will  be  used  in  this  paper  to  determine 
the  Cr . 
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3.  PROGRAM  OVERVIEW 

Estimation  of  the  parameters  for  use  in  the  Prohl 
sequence  calculations  is  a  choice  made  by  the  author.  The 
values  of  m,  A,  I  and  J  for  each  section  could  be  obtained 
from  scantling  lists  and  a  weight  curve  but  their  assembly 
would  be  time  consuming  and  defeat  one  of  the  purposes  of 
this  work,  that  is  rapid  vibration  information  with  a  minimum 
of  data  required  as  input.  Therefore  the  subroutine  HULL  for 
estimating  these  quantities  has  been  incorporated  into  the 
program.  Input  required  is  station  spacing,  number  of 
stations,  displacement,  thickness  of  side  shell,  main  deck 
and  bottom  plating,  offset  data  and  deck  edge  heights  for 
each  station. 

Sectional  mass  is  computed  by  first  obtaining  the 
total  sectional  area  from  the  offset  data  for  each  station. 

A  partial  sectional  area  is  computed  by  the  use  of  the  five- 
eight  rule  [14]  and  the  offset  data.  An  additional  area  is 
computed  which  is  equal  to  the  deck  edge  height  minus  the 
highest  waterline  for  which  offset  data  is  input  times  the 
offset  at  that  waterline.  The  sectional  area  equals  the 
sum.  A  uniform  permeability  is  a  reasonable  assumption  for 
naval  vessels  so  that  the  mass  is  proportional  to  the  area 
at  each  station.  The  sum  of  the  sectional  areas  is  divided 
into  the  displacement  to  get  an  area  density.  The  density 
is  then  multiplied  by  each  sectional  area  to  obtain  a 
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sectional  mass.  No  attempt  is  made  to  account  for  super¬ 
structure  concentration  of  mass  or  its  effect  on  stiffness. 

A  two  beam  ship  model  would  be  the  next  step  in  sophistication 
The  mass  moment  of  inertia,  J,  about  the  neutral  axis 
is  determined  by  first  assuming  that  the  beam  cross  section 
has  a  uniform  distribution  of  mass.  Naval  vessels  are 
characteristicly  longitudinally  framed.  In  the  midship 
section  the  hull  strength  is  provided  by  the  main  deck, 
side  shell  and  bottom  plating  and  their  associated  longi¬ 
tudinals.  A  thirty  inch  spacing  is  not  uncommon  and  the 
stiffener  area  to  plate  area  ratio  is  about  .3.  So,  the 
effect  of  longitudinal  is  taken  into  account  by  increasing 
the  plating  thickness  by  30%.  The  neutral  axis  is  determined 
by  a  calculation  of  the  first  moment  of  the  hull  steel.  The 
sectional  area  density  is  then  used  along  with  the  offset 
data  in  a  second  moment  calculation  about  the  neutral  axis. 
Each  hull  cross  section  is  divided  by  waterlines  with  four 
foot  spacing.  The  area  from  the  baseline  to  the  first  water¬ 
line  is  assumed  to  be  bounded  by  a  parabola.  All  subsequent 
areas  are  assumed  to  be  trapezoids.  The  distance  from  the 
neutral  axis  to  the  centroid  of  each  incremental  area  is 
computed,  squared  and  multiplied  by  the  area.  The  second 
moment  of  the  area  about  its  own  centroid  is  computed.  The 
sum  of  all  these  terms  equals  the  sectional  second  moment. 

This  is  multiplied  by  the  sectional  area  density  to  obtain  J. 


■fTUgjg,"1  M  *E 


The  moment  of  inertia,  I,  is  obtained  by  a  second 
moment  calculation  of  the  hull  steel  around  the  neutral 
axis. 

v 

The  effective  shear  area  KA  in  a  hull  shape  is  very 
difficult  to  determine.  Traditional  practice  has  been  to 
use  only  the  vertical  side  shell  and  any  continuous  longi¬ 
tudinal  bulkheads  in  the  calculation.  Analytically  obtained 
values  for  K  are  for  simple  symmetric  shapes  only.  [15] 
Calculations  by  the  author  for  box  type  sections  with 
various  side,  deck  and  bottom  thickness  and  beam/draft  ratios 
shows  that  an  approximate  value  for  the  effective  shear  area 
can  be  obtained  from  the  following 

(KA) i  =  . 87 ( 2tD)  =  1 . 7  4  tD  (1) 

where  t  is  the  sideshell  thickness  and  D  is  the  depth.  Then 
for  steel 


KAGi  -  1. 74tDE/2 (1+u) 


or 


KAGi  »  2.35  x  108  tD 


(2) 


for  t  in  inches  and  D  in  feet. 


Subroutine  NMODES  first  computes  the  virtual  mass 
of  each  station  which  varies  with  the  mode  of  vibration  for 
which  a  solution  is  sought.  Burrill's  formula  [14]  is  used 
to  estimate  the  frequency  of  two  mode  vibrations  and  the  J 
factor  for  that  mode  is  formed  and  virtual  masses  computed. 
Equations  (1)  through  (4)  of  Section  B  are  solved  for  each 
station  given  a  value  of  oj  in  the  following  way:  the  shear, 
bending  moment,  displacement  and  slope  at  any  station  i  are 
assumed  to  be  of  the  form 


(3) 


V.  =  a .  +  b. 9 
i  i  is 


M.  *  -c.  +  d. 9_ 
i  i  is 


y.  =  e.  +  1^9, 


0.  =  g.  +  h. 0_ 
i  i  is 


where  9  is  an  unknown  initial  slope  at  the  stern.  Then  (1) 

5 

through  (4)  of  Section  B  are  used  to  determine  the  conditions 
at  station  i+1 


ai+l  *  ai  +  mt“  ei 


bi+l  ”  bi  +  mi“  fi 
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'i+1 


g.  -!ii 

yl  EI , 


ai+1L 


i+i 


hi+l 


bi+1L 

hi  EI, , ,  2EI , 


•i+1 


i+1 


"ci+l  “  "c  +  a^iL  *  “  5, 


i+1 


i+1 


'i+1 


di+l 


di  +  bi>iL  -  Ji+i  "2  hi+i 


(4) 


•i+1 


fi+l 


ciL 

ei  +  giL  "  2Ei: 


ai+lL' 


dLL 

fi  +  hiL  +  HIT 


i+1 

2 


6EI 


i+1 
3 


1+1 


bi+lL 
6EIi+l 


ai+lL 

KGAi+i 


bi+lL 

KGAi+i 


From  the  treatment  in  Section  B 


0  = 
s  d 

N 


(5) 


and 


N 


Vn  .  Q 


(6) 


at  a  natural  frequency.  The  subroutine  performs  the  sequence 
calculation  until  opposite  signs  are  obtained  for  the  quantity 
on  the  left  hand  side  of  (6)  for  two  values  of  a.  A  natural 
frequency  is  then  known  to  lie  between  these  two  values. 
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The  subroutine  is  designed  to  determine  the  natural  fre¬ 
quencies  to  an  accuracy  of  ±.04  radian/sec.  After  a  natural 
frequency  is  obtained  the  mode  shape  is  determined  from 


y  »  e  +  f 
1  i  1  dN 


(7) 


and  normalized  by  the  method  discussed  in  Section  B 


N 

Z 

i=l 


m. 


i 


1 


Data  obtained  from  the  Andersson  reference  [2]  for  a 
60,000  ton  tanker  was  input  into  NMODES.  The  resulting 
natural  frequencies  obtained  are  compared  with  the  measured 
frequencies  in  radians/sec. 


NMODES 

[2] 

1 

4.92 

4.89 

+  .6 

2 

10.63 

10.43 

+1.9 

3 

16.68 

16.41 

+1.6 

4 

23.20 

22.91 

+1.3 

5 

29.26 

29.17 

+  .3 

6 

34.88 

35.03 

-  .4 

7 

39.77 

40.33 

-1.4 

percent 

difference  is  quite 

small . 
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The  input  subprogram  reads  data  which  is  punched 


according  to  the  format  listed  after  each  read  statement. 

The  data  must  be  placed  in  the  proper  fields,  decimal  points 
must  be  punched  and  all  integers  must  be  right  justified. 

STA  is  an  integer  variable  and  is  the  number  of  stations 
for  which  offset  data  is  given.  STASPA  is  the  station 
spacing  in  feet,  DISPL  is  the  displacement  of  the  vessel  in 
long  tons.  The  program  later  changes  this  value  and  does 
all  calculations  in  slugs,  BEEM  is  the  beam  of  the  vessel 
in  feet.  (Another  variable  with  the  correct  spelling  is  used 
in  the  program.)  Tl,  T2  and  T3  are  the  main  deck,  side  shell 
and  bottom  plating  thicknesses  respectively  of  the  midship 
section  in  inches.  DECXE  is  the  deck  edge  height  from  bow 
to  stern  in  feet.  There  must  be  one  entry  per  station.  DATA 
is  offset  data,  one  card  per  station.  All  entries  are  in 
decimal  feet.  The  first  field  is  a  real  integer  (decimal 
point  must  be  included)  which  is  closest  to  the  draft  at 
that  station  divided  by  4.  The  remaining  twelve  fields  are 
offset  data  starting  from  the  bottom  and  proceeding  upward 
at  four  foot  intervals.  The  first  field  must  have  a  real 
integer  less  than  or  equal  to  twelve  and  only  12  waterlines 
may  be  input,  unused  fields  may  be  punched  with  zeros. 
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The  output  will  list  the  computed  station  mass  in 

2 

slugs,  the  mass  moment  of  inertia  J  m  slug  ft  ,  the  second 

4 

moment  of  area  I  in  ft  ,  shear  stiffness  KAG  in  lbf  and 
the  2-D  added  mass  for  each  station  in  slugs.  The  program 
will  list  the  first  seven  natural  frequencies  and  the 
normalized  made  shapes.  The  program  will  also  list  the 
response  amplitude  operator  over  the  frequency  range  of 
interest  for  speeds  of  0,  10,  20  and  30  knots.  It  will  also 
list  the  response  of  the  vessel  to  a  Pierson-Moskowitz 
fully  developed  sea  with  a  wind  speed  of  50  knots.  Response 
amplitude  is  in  inches  measured  at  the  bow. 

It  is  obvious  from  a  look  at  the  Pierson-Moskowitz 
formula  that  precious  little  energy  exists  at  frequencies 
above  the  first  mode  of  vibration,  springing.  It  is  reason¬ 
able  then  to  assume  that  wave  excitation  is  significant  only 
for  large,  long  and  shallow  vessels.  The  large  masses  and 
the  low  stiffnesses  (due  to  reduced  depth)  are  responsible 
for  a  relatively  low  first  natural  frequency,  one  low  enough 
that  sufficient  energy  is  available  for  excitation.  Too,  as 
wave  frequency  goes  up  the  effect  on  a  vessel  of  constant 
draft  drops  off  logarithmicly  as  has  been  shown.  The  only 
times  then  that  the  higher  modes  will  be  excited  are  when 
slamming  occurs  or  under  some  other  impulse  type  loading. 

The  value  of  the  program  lies  in  the  beam  parameters  obtained 
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from  the  HULL  subroutine,  the  natural  frequencies  and  mode 
shapes  obtained  from  NMODES  and  the  forcing  data  which  can 

be  retrieved  from  FORCE.  The  node  location  in  the  natural 

•  v 

modes  is  important  as  possible  locations  for  placement  of 
machinery  operating  at  or  near  that  natural  frequency  since 
excitation  at  a  node  will  not  excite  vibration. 

Numerous  comment  cards  have  been  inserted  in  the 
program  and  can  be  used  to  understand  the  many  small  pro¬ 
cedures,  which  are  a  natural  part  of  any  project  of  this 
dimension.  The  text  of  this  thesis  serves  only  to  discuss 
the  major  technical  areas  and  not  the  interesting  but 
obviously  minor  structural  details  of  software  construction. 
The  program  listing  which  follows  is  divided  into  the  sub¬ 
routines.  A  sample  input  record  and  output  record  are  shown 
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0 . 2423 
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0.2742 

06 

0.7262 

93 
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0. 2o9a 
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0 . 5962 

08 

5 

0.  2  322 

0  6 

0.6262 

03 

6 

0.290a 

06 

0 .6432 
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0.29  42 

2  6 

C . 6532 

33 

8 

0. 2962 

0  6 

0.6533 

08 

9 

0.2972 

06 

0 . 660  3 
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10 

0. 2972 

9  6 

0  .  6  6  0  S 

03 

1 1 

0.2972 

06 

>• 

0.6602 

03 

12 

C.2972 

Oo 

0.6603 

ca 

13 

0. 2972 

96 

3. 6602 

03 

14 

0.2962 

06 

3.6532 

33 

15 

0.2942 

0  6 

0.653a 

0  3 

16 

9.2902 

9  o 

0 . 6  4  3  2 

36 

17 

0.  2322 

06 

0 . 6  2  n  2 

0  3 

18 

0.269  2 

0  6 

3.  5962 

03 

19 

0. 24cj2 

0  6 

0. 5512 
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2  0 

9.2422 
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0 .6422 

36 
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Ai)  0  2  AS 

0.2722 
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3 . 2b  72 

1  1 

3,150.1.  j 6 
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1  1 
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3.2262 

1  1 

2.333a  0  i 

3.2672 
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3.  27  32 
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3 .  zz  b  a 
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0  5 

3.2262 

1 1 

3.  4052  06 
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05 

0.2262 

1  1 

0.4052  06 

0.2792 

35 

3.2263 

1 1 

0 .  40  52  0  6 

0.2792 

G5 

9.2262 

1  1 

3  .  4  352  3 o 

0.2792 

05 

0.2262 

1 1 

0.4052  36 

3.2762 

0  5 

3 .2262 

1  1 

3 .  40  3a  3  6 

0. 2772 

05 

3.22oa 

1 1 
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0.2752 

35 

0  .  2  2  o  2 
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1  1 

3 . 36  62  0  : 

C.  2562 

*  £ 

'j 

3 . 2  2  o  2 

1  i 

•  j  3  j  i  *.  o 

3.2392 

3  5 

■3  .2262 

1  1 

j  •  ^ 3 o m 

3.2642 

05 

0.2472 

1 1 

3.  2242  3o 

Sample  Output 


106 


i mz.' . - - - 


NATURAL  FREQUENCY  NU.135E  1  IS  2.407  HAD/ SEC 

0.578 jE-03  0.44162-03  3.30392-03  0.18321-03  5.39252-'4 

-0.77213-04  -  0.  176  IE-03  -0.25271-03  -0.30303-03  -0.32452-03 
-0 .3 15 9 £-03  -C.277SE-03  -0.21242-03  -0.  12302-03  -3.  14473-'* 4 
0.  10312-03  0.23922-03  0. 37332-03  0.50602-03  5,o394>33 

NATURAL  F3EQ0ZNCT  NU3BZ5  2  IS  6.040  EAL/SZC 

0.5409E-03  0.32442-03  0. 10622-03  -0.94752-04  -0.25342-03 

-0  .34o  5ii-C 3  -0.36302-03  -0.  30402-03  -0. 18351-03  -0.26332-04 
0.13592-03  0.27073-03  3.35052-03  C.35792-G3  0.23762-53 

0.14862-03  — J. 40602-04  -0.25552-03  -0.4709--J3  -0.66992-03 

NAT  DEAL  FREQUENCY  N0I132R  3  IS  10.611  EAD/S2C 

0. 49662-03  0.21502-03  -0.63352-04  -0.23C6E-03  -:. 27482-03 

-0.31832-03  -0.13752-03  0.98632-04  0.30222-03  0.39612-03 

0.35112-03  J. 17872-03  -0. 55202-04  -0.2o342-03  -0 . 3 667 2-0 3 
-0 .32253-03  -O.1 3812-03  0.  1351E-03  C.  42532-03  0.57182-03 

NATURAL  FREQUENCY  NU.1B2R  4  IS  15.573  EAD/SEC 

C. 45472-03  0.11382-03  -0.20532-03  -0. 37652-03  -0.31532-03 

-0.58122-04  0. 23992-03  0.39932-03  0.32322-03  0.55562-04 

-0 .24543-03  -o. 40062-03  -C. 319 12-03  -0.51192-04  5.^4032-03 

0.37732-03  0.27122-03  -0.31992-04  -0.33962-03  -0.o6082-03 

NATURAL  FREQUENCY  NU.1322  5  IS  2C.  577  EAC/SZC 

0.41912-03  0  .23392-04  -0.3101E-03  -0.37192-03  -0.  1  1  162-03 

0. 2551E-03  9.40972-03  0.2133E-03  -0.16622-03  -0. 39852-03 

-0.28092-03  0.8392  E-04  0.37762-03  0.  3450E-03  0.16002-04 

-0.31312-03  -0.35562-03  -0.61532-04  0.35852-03  >.o445Z-5'3 

NATURAL  FREQUENCY  NU .1 3 2 E  6  IS  25.  436  5AD/S2C 

0. 38742-03  -0.57652-04  -0. 37482-03  -C. 27732-  53  0.l4ei--:3 

0.41702-03  0  .2  09  32-03  -0. 24342-03  -0.40232-03  -  0. 8450E-'4 

0.33402-03  0.35452-03  -0,48452-04  -0.39672-03  -0.27962-03 

0.15712-03  0.39202-03  0.  1463  2-63  -0  .33322-0  3  -0.  o  5  77  :-,;3 

NATURAL  FREQUENCY  NUN3ER  7  IS  29.985  EAD/SEC 

0.35452 -0  3  -0.12602-03  -3.39352-u3  -0.  12212-03  0. 3  ■♦89  2-03 

0.32762-03  -0.17332-03  -  0.40792-03  -0.19722-04  0.40  .02-03 

0.21062-03  -0.29892-03  -0. 35052-03  0.13552-53  0.421 1 2-53 

0.806o2-04  -C. 36122-53  -0.241  32-03  0.28882-03  0.  59592—53 
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C<  O  v  > 


V  =  0  7*10 


??.  ZQ  UZSC7 

?.  AC  5 

ZSPQUSZ 

5.  AO 

?.  Z  3?  OH  S  2 

1.2596 

4 . 2  o  5 

0.4529 

2  4,64 

1.  232 

0.277  0 

5.  9o  8 

0.7934 

26.52 

2.10  a 

0.2945 

6  .  )  1  6 

1.213 

23.30 

2.  3c 0 

0.3119 

10.32 

1.079 

25. 03 

3 .  c  3  1 

0.3293 

12.77 

2.  120 

28.  79 

4.142 

0 . 3  46  8 

15.27 

2.533 

27.34 

4.  259 

0.3642 

17.72 

2 .  do  1 

24.49 

3.991 

0.3817 

19.99 

3.096 

2C.1  1 

3.346 

0.3991 

21.35 

3.227 

14.32 

2.  336 

0 . 4  16  5 

23. 46 

3.253 

8.22  0 

1. 356 

0.4340 

2  4 . 4  0 

3.176 

3.562 

1. 333 

0. 4514 

24.66 

3.003 

17.34 

2.721 

C .4669 

24.  17 

2.  746 

25.26 

4.283 

9.4363 

2. «  o d 

2.427 

39.49 

5.766 

0.5C37 

20.3  1 

2.0  59 

4  9.  9  5 

7.  001 

0.5212 

16.04 

1 . 664 

53.59 

7 . 3b5 

0 . 538  6 

14.70 

1 .265 

64.37 

.9 . 2  o  o 

0.5560 

10.93 

0.6b2d 

66.3  3 

3.  145 

0.5735 

7.  134 

0.5359 

63 . 5  0 

7.438 

0 .590  9 

3.  43C 

0.2410 

56.58 

6.331 

0.6084 

0.5410 

0.3  5602-C  1 

44.57 

4.763 

0.6253 

2.  690 

0.1659 

28.67 

- .  925 

0.6432 

4.  545 

0.  2631 

10.31 

1.054 

0.6607 

5.501 

0.2993 

11.73 

1 . 0"5  3 

0.6781 

5.  552 

0.2341 

29.24 

2.  605 

0.6956 

4.605 

0.2316 

43.73 

3.726 

0.7130 

3.  466 

0.1575 

52.54 

4. 233 

0.73C4 

1.  3C9 

0.7762E-01 

54.  10 

4.230 

0.7479 

0 .1637 

0.6340S-U2 

47.9) 

3.590 

0.765  3 

1.  2a0 

0.491  U-ul 

34.00 

2.492 

0.78 26 

2.207 

0.30162-u  1 

16.  78 

1.  156 

0.8002 

2. 545 

0 . 3764  c-0 1 

3.525 

) .  5  a  5  3 

0 . 6 17o 

2.  307 

0.75402-01 

25.  43 

1.621 

0.6351 

1.622 

0.50362-01 

39.  06 

2.395 

0.8525 

0. 7051 

0.  20322-0 1 

44.73 

4,  •  o  3 

0 . 869  9 

J.2318 

C.  6  51  4Z-02 

41.  10 

?  i  *  6 
2*  • 

0. 8374 

0.3*9) 

0.2408 2-0 1 

26.74 

1.575 

0.9046 

1.  215 

0.31001-01 

11.57 

0.6  115 

0.9223 

1.142 

0. 275 12-01 

13.62 

0. o945 

0. 9397 

) . 7693 

0. 17902-01 

29. bo 

1  .  4  o  1 

0.9571 

4.2937 

0 . 6  5  2  3  2-  0  2 

3  o .  3  2 

1.  a47 

0.9746 

0.3  054 

0 . 6  49  22-02 

37.70 

1.  7  34 

0.3920 

0 • 5443 

0.11)72-01 

26.50 

1.179 

1  .CC9 

0.54*3 

1  07  22-01 

10.05 

).  4326 

1.027 

0.3  514 

0.65662-02 

1 5,  o3 

0.3541 

1.044 

0.  1576 

0 . 2625&-02 

30.52 

1.23* 

1.062 

9.2312 

0.  54092-02 

36.33 

1.426 

1.  37  9 

4 . J  596 

). 59412-02 

30.6  0 

1 .153 

1.097 

0.2442 

0. 36772-02 

15.54 

0.5731 

1.  114 

0,52322-01 

0.79332-06 

10.26 

0. 3632 
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ft"' *3** 


V  =  2 )  7  =  30 


EAO 

E35?C:iSL 

a  AC 

5.  IS  PC  NS  1 

52.31 

1.  0  J9 

79.66 

0.5709 

53.28 

2.253 

69.32 

1.  472 

63.44 

3.619 

96.67 

2.574 

67.29 

5.  570 

105.6 

5.005 

69.27 

7.213 

109.9 

7.301 

66.33 

3.  450 

110.5 

9.467 

65.42 

9.0  25 

106.6 

11.06 

58.60 

9.767 

97.41 

1 1.68 

48.  13 

7.  599 

32.57 

11.02 

34.11 

5.563 

o  1.96 

3.951 

17.89 

2.971 

36.54 

5.550 

13.36 

2.231 

15.  57 

2.  478 

32.81 

5.454 

37.64 

6.  148 

56.21 

9. 237 

7  2.  96 

12.10 

79.11 

12.73 

103. 3 

13.07 

99.  14 

15.67 

139.7 

23.29 

114.  1 

17.57 

163.7 

27.  13 

121.8 

18.25 

177.2 

29.04 

120.7 

17.54 

177.  4 

26. 65 

1C9.7 

15.46 

162.9 

25.8  5 

89.03 

12.  14 

133.5 

20.76 

59.83 

7.889 

91.01 

13.35 

25.40 

3.236 

40.07 

5.955 

13.50 

2.  277 

25.38 

3.680 

52.96 

6.297 

76.68 

1C. 84 

82.65 

9.492 

121.6 

16.73 

101.6 

11.27 

150.7 

20.19 

106.4 

11  .40 

158.7 

20.69 

95.66 

9.  903 

143.4 

13.19 

70.56 

7.06C 

106.5 

13.  13 

35.54 

3.438 

54.31 

6.515 

15.42 

1. 442 

22.39 

2.612 

48.53 

4.396 

71.73 

'6.137 

76.50 

6. 700 

H3.9 

12.57 

68.86 

7.  535 

132.3 

1  *  .  2  7 

32.41 

6.768 

123.7 

12.92 

58.  38 

4.  646 

33.  32 

3.943 

24.12 

1.361 

36. 63 

3.  bz3 

26.24 

1.963 

3  6.86 

3.742 

53.55 

4.  251 

37.  44 

8.  193 

77.39 

5.453 

116.0 

10.59 

75.64 

5.  175 

113.6 

10.10 

53.54 

3.559 

a  0.53 

6.  987 

20.54 

1.327 

31.03 

2.523 

31.04 

1 . 9  50 

46.  40 

3.324 

6C.94 

3.7  24 

91.37 

7.  346 

72.93 

4.  336 

109,5 

6.533 

61.48 

3.558 

92.  36 

7.071 

31.29 

1.764 

47.05 

3.517 

20.6  2 

1.  132 

30.95 

2.259 

Sample 
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